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A hallmark feature of topologically ordered states of matter is the dependence of ground state 
degeneracy (GSD) on the topology of the manifold determined by the global shape of the system. 
Although the topology of a physical system is practically hard to manipulate, recently it was shown 
that in certain topologically ordered phases, topological defects can introduce extra topological 
GSD. Here the topological defects can be viewed as effectively changing the topology of the physical 
system. Previous studies have been focusing on two spatial dimensions with point-like topological 
defects. In three dimensions, line-like topological defects can appear. They are closed loops in the 
bulk that can be linked and knotted, effectively leading to complex three dimensional manifolds 
in certain topologically ordered states. This paper studies the properties of such line-defects in 
a particular context: the lattice dislocations. We give an analytical construction, together with 
support from exact numerical calculations, for the dependence of the GSD on dislocations of certain 
doubled versions of the exactly solvable Kitaev's toric code models in both two and three dimensions. 
We find that the GSD of the 3d model depends only on the total number of dislocation loops, no 
matter how they are linked or knotted. The results are extended to Z n generalizations of the 
model. Additionally, we consider the phases in which the crystalline orders are destroyed through 
proliferation of double dislocations. The resulting phases are shown to host topological orders 
described by non-Abelian gauge theories. 
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I. INTRODUCTION 

The introduction of topological order as the underlying 
principle of the fractional quantum Hall effect (FQH) 1 
was central to the development of this new paradigm in 
modern condensed matter physics, which parallels the 
ubiquitous symmetry breaking mechanism. Topologi- 
cally ordered phases, by definition, are quantum states of 
matter having long-range quantum entanglement, which 
can be characterized via the topological entanglment 
entropy 23 . Topologically ordered phases can exist in 
both two and three spatial dimensions. 

In topologically ordered phases, non-trivial excitations 
such as fractional charges may arise and realize any- 
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onic statistics in two spatial dimensions^. Tunneling of 
virtual pairs of such fractional charges across the sys- 
tem can lead to a degenerate ground state, which im- 
plies non-trivial ground state degeneracy (GSD) that de- 
pends on the topology of the underlying physical sys- 
teirP, which is another hallmark of topological ordeiR A 
system with excitations is typically also described by de- 
generate states, and spatial manipulations of excitations 
correspond to unitary transformations (Abelian or non- 
Abelian) of the degenerate states,potentially realizing 
quantum computational operations'"^ Recent years saw 
experimental implementations of anyonic excitations and 
other ingredients of quantum computation in photonic, 
superconducting, ultracold atoms, and most promisingly 
in non- Abelian FQH systems [for a review, see Refs l9ll(?] . 

In particular, the specific dependence of the ground 
state degeneracy on the real space topology of the phys- 
ical model is a crucial signature of topological order, 
but is hard to probe directly^. Recently, a novel pro- 
posal to circumvent this constraint was made 12 for a 
FQH system on a (2d) lattice. That proposal is based 
on the mapping of a flat band with Chern number N 
onto N independent quantum Hall systems ("layers"), 
while these layer degrees of freedom are connected by 
the two-dimensional lattice translations. Lattice disloca- 
tions (translational topological defects) are then shown 
to act as ribbons connecting the layers, effectively setting 
the quantum Hall system on a higher genus surface, al- 
though the physical system shape is unchanged. Further, 
the dislocations themselves carry a non-trivial quantum 
dimen sion e ven in an Abelian FQH system^. i n earlier 
studie^ES it was also shown how lattice dislocations in 
the 2d Toric Code model can change the anyonic exci- 
tations of the model. This series of works demonstrate 
qualitatively new features of topological defects (point- 
defects) in topologically ordered phases in two spatial 
dimensions; namely these defects can carry extra GSD. 

However, in three spatial dimensions, topological de- 
fects can be either point-like or line-like. In three spatial 
dimension, a new ingredient in the problem is that the 
line-like defects can be linked or knotted. For example, 
crystalline dislocations in a three-dimensional lattice are 
line-like defects which form closed loops. Therefore a 
dislocation itself can have a complicated topology if it is 
knotted or linked with other dislocation loops. 

Natural questions emerge: Are there qualitatively new 
features associated with line- like topological defects (e.g., 
dislocations) in 3d topologically ordered phases? Can the 
line-like defects carry extra GSD? In addition, because 
the GSD is a topological property of the system, GSD 
can only be a topological invariant of the links and/or 
knots of the dislocation loops in 3d. Is the GSD a trivial 
or non-trivial topological invariant of the links/knots? 

It is interesting to speculate what would happen if 
the GSD was indeed a non-trivial topological invari- 
ant of the links/knots. In this case, the non-trivial 
topological invariant is distinct from the conventional 
topological invariants of links/knots (such as the Jones 



polynomia l 15 * 16 *) : The conventional topological invari- 
ants of links/knots are mappings from links/knots to a 
complex number. But here, the possible non-trivial topo- 
logical invariant of the links/knots is a mapping from 
links/knots to a Hilbert space associated with the GSD. 
We note that this new type of topological invariants of 
links/knots has attracted some interest in mathematics^ 
and physicsIS 

At this point, it is worth mentioning that qualita- 
tively new features of topological defects in topological 
phases without long-range quantum entanglement (and 
thus not topologically ordered) were shown quite some 
time ago. Examples include Majorana fermions bound 
to vortices of p-Mp topological superconductors^, and 
helical modes bound to dislocations in 3d topological 
insulators^. Even general understanding for such phe- 
nomena is available in the context of non-interacting 
fermion systems 21 . However, our understanding of topo- 
logical defects in topologically ordered phases is probably 
far from complete. 

Motivated by these questions, in this paper we attempt 
to obtain a general understanding, to a certain level, of 
the extra GSD induced by dislocation loops in certain 
topologically ordered phases in 3d. We mainly focus on 
a 3d topologically ordered phase described by Z 2 x Z 2 
gauge group, in which the lattice translation symmetry 
interchanges the two Z 2 gauge sectors. Generalization to 
other topologically ordered phases with Z n x Z n gauge 
groups is also presented. To study the effects of disloca- 
tion defects in these 3d phases, we use exactly solvable 
models as a powerful tool. 

In the past, important insights about topological or- 
der came from class es of exactly solvable models in 
twcEEUMM] and thretPMHH! dimensions. The paradig- 
matic exactly solvable model is the two-dimensional 
"toric code"^ (2dTC) defined through interactions of 
spin-| sites on a square lattice. This model was used to 
originally introduce the concept of topological quantum 
computation using anyons^. Signatures of its Abelian 
anyonic excitations have been observed and further ma- 
nipulation schemes are developed**"^. Parallel to such 
efforts, the toric code has been generalized to three 
dimensions^ 1 (3dTC). Both 2dTC and 3dTC are de- 
scribed at low energies by a Z 2 gauge theory, describ- 
ing therefore Z 2 topological order. However, the 3dTC 
is predicted to stay topologically ordered even at finite 
temperatures due to its dimensionality^. The 3dTC is 
still a relatively simple model, and therefore we can study 
in detail its properties on the three-dimensional lattice. 

We calculate the effect of arbitrary dislocation loops 
in an exactly solvable three-dimensional doubled Kitaev 
toric code^. The model system consists of two copies 
("flavors") of the 3dTC, displaced by half of lattice con- 
stant along one direction (say y), which allows dislocation 
lines with the minimal Burgers vector b — y/2 to connect 
the two flavors (sublattices) and therefore induce non- 
trivial topologies of the total (three-dimensional) man- 
ifold on which the model exists. We find that our con- 
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structive analytical results, which are corroborated by ex- 
act numerical calculations, show that this Abelian model 
is not sensitive to the knotting and linking of dislocation 
lines, but can only distinguish the total number of dislo- 
cation loops in the crystal lattice. The topological ground 
state degeneracy in a lattice with periodic boundary con- 
ditions scales simply with the number of dislocation loops 
(k ^ 1) as \GSD\ = 2 6 • 2 k ~ l (\GSD\ = 2 6 when disloca- 
tions are absent). We generalize this result to the dou- 
bled Z„ toric code in three dimensional periodic lattice 
in the presence of dislocations, where the ground state 
degeneracy is \GSD\ = n e ■ n k ~ 1 , again independent of 
linking and knotting of the dislocation loops. This im- 
plies a zero-temperature entropy contribution which is 
extensive in the number of dislocation loops, and we will 
show it seems to persist in a realistic system with open 
boundaries (i.e. without periodicity). 

We find somewhat similar results in the two- 
dimensional doubled Kitaev toric code, where a dislo- 
cation loop takes the form of a pair of pointlike edge 
dislocations (with Burgers vectors b and —b). The GSD 
of the lattice with periodic boundary conditions (equals 
2 4 in absence of dislocations) is also not influenced by a 
single dislocation pair, but scales with number of pairs 
k > 1 as \GSD\ = 2 4 • 4 k -\ Further, \GSD\ = n 4 • n 2k - 2 
for doubled Z„ toric code in periodic system in two di- 
mensions. 

A further fundamental question arises: What is the 
state of matter obtained by making the dislocations dy- 
namical excitations? In 2d, we consider a state where 
"double dislocations" (|6| = 1, therefore not mixing fla- 
vors) proliferate and melt the lattice. Then, a tunneling 
event of a b — y/2 dislocation pair introduces a Z2 twist 
across which the flavors are exchanged. We show that 
the resulting theory is described by a non- Abelian group 
G n = (Z„ x Z n ) x Z2 by analyzing the ground state de- 
generacy on the torus. We further derive the degeneracy 
on arbitrary genus oriented surfaces, which reveals the 
quasiparticle content of the theory. In the case n = 2, i.e. 
the melting of simplest doubled toric code, both the de- 
generacy and the quasiparticle content explicitly match 
the properties of topological states with (non-Abelian) 
dihedral group D4, as is already indicated by the isomor- 
phism D4 ~ (Z2 x Z2) x Z2. In 3d, we focus on the melted 
state on a three-torus, i.e. a three-dimensional melted 
lattice with periodic boundary conditions. Also, we only 
consider the original toric code, i.e. n — 2, and find 
the ground state degeneracy of the melted state, which 
exactly matches the degeneracy of the D4 topologically 
ordered state on the three-torus. Therefore, there is very 
strong evidence that the condensation of double dislo- 
cations, which carry no flux nor charge in the Abelian 
Z„ x Z„ theory, gives rise to a non-Abelian gauge theory 
in both two and three dimensions. The defect conden- 
sation promotes a global symmetry (the Z2 flavor) to 
a gauge symmetry, and can be seen as a physical way 
to connect phases that exhibit different topological or- 
der and different global symmetry, i.e. connect different 



symmetry enriched topological phases. 

This paper is organized as follows. The majority of 
the discussion is focused on the doubled Kitaev toric 
codes, and we leave the detailed analysis of the gener- 
alized doubled Z„ toric code for Appendix [A"| We start 
with the complete analysis of the doubled 2dTC model 
in Section [TTJ The model description including disloca- 
tions, and proof of the ground state degeneracy are pre- 
sented in detail, and also serve as a warm-up for the 
3d case. The exact numerical results for the 2d model 
are presented jointly with the 3d case in Section |IV A| 
In Section imi we introduce the three-dimensional lattice 
model and construction of dislocation loops. Some de- 
tails of the construction are relegated to Appendix[E] We 
present the behavior of stabilizer constraints, the string 
and membrane operators, and the exact numerical results 
for GSD in Sections |IVA[ |IVB| Section |IVC| introduces 
the rules for string operators we found by analyzing the 
lattice. This facilitates the presentation of our main ana- 
lytical results for arbitrarily linked dislocations with, and 



without, knotting in Sections IV D and IV E[ respectively. 
The detailed analysis from the dual perspective of mem- 
branes is presented in Appendix [B] Section [V] presents 
the analysis of the theory of melted lattice, with emerg- 
ing non-Abelian gauge theory. We close with a discussion 
of Z„ models (which is studied in detail in Appendix |A|, 
the influence of system boundaries, as well as some im- 
plications of our study. 



II. THE TWO-DIMENSIONAL DOUBLED 
TORIC CODE 

The 2dTC is a very well knowrPl model that has Z2 
topological order. We describe our doubled variant in 
detail below. 

The square lattice of one 2dTC copy has lattice con- 
stant a = 1 and spin— |'s on lattice links. The second 
copy of the 2dTC is displaced by half the lattice constant 
along y direction, this direction being vertical in Fig. [T] 
The two copies are called "flavors", having flavor index 
/ = a, (3. When we discuss one copy of the model, i.e. 
one fixed flavor, it is understood that we consider the 
copy of the square lattice assigned to that copy, having 
a = 1. Star operators (circles of the same color in Fig. [IJ 
are positioned at lattice sites and act only on the four 
neighboring spins (Fig. [T]) : 
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J s+x/2 (J s+y/2 s-x/2 c 



(1) 



where o~ l R is the i-th Pauli matrix acting on the spin at po- 
sition R, while s is a lattice site position in the / flavored 
lattice. Each square plaquette of the lattice centered on 
p = s + (x/2 + y/2)/2 for some s, carries a plaquette 
operator B acting on 4 neighboring spins: 

Bl 



7 p+x /2°p+y /2 a p-x /2 a p-y /2 



(2) 



The centers of plaquettes belonging to flavor a(/3) are 
easily recognized as positioned on spins of flavor /3(a), 
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FIG. 1: The 2d doubled toric code with a dislocation pair. Red/yellow dots are spins of flavor a//3, two copies displaced 
by ay/2. Green/blue circles mark centers of a//3 star operators. Removing all sites within the shaded area (the Volterra, 
u aj3" , line) leaves two dislocations (positioned at blue shaded sites). A s is a trivially locally repaired star operator (acting 
with o x on neighboring spins labeled by wavy lines). Near the dislocation site, non-trivial local surgery reinstates commuting 
stabilizers: Star marked by black triangle now acts on three light blue diamond spins only; Two plaquettes marked by black 
triangles are replaced by the merged 9-spin stabilizer (acting on violet square spins). If the dislocation was sitting at a star 
operator site (example not shown), the two new stabilizers would look the same, except that the assignment star/plaquette 
would be exchanged. Dashed line follows a lattice string loop P, giving operator £2^ acting through a z on crossed-out spins. P 
is constructed locally step-by-step, always hopping between two star lattice sites and acting on the spin site they share. Notice 
the flavor change (flavor is also marked by slant of cross-out lines) when crossing the a/3 line, leading to the obligatory even 
winding around the dislocation. Multiplying E2d with the plaquette B p locally changes the shape of the loop P, but not the 
(possible) independence of string loop E2d from the set of plaquette stabilizers. The T,2d is actually trivial, because the loop is 
contractible. One must explicitly use the repaired plaquettes near the dislocation site to show this. 



when there is no star operator (circle) on that spin al- 
ready, see Fig. [T] The displacement of two copies by 
y/2 (vertically) creates alternating columns of star and 
plaquette operators in the total 2d lattice. 
The Hamiltonian 

h = - E fe^+E**) (3) 

f£{a,/3} \ s V ) 

is formed by the star and plaquette operators ("stabiliz- 
ers"), defined on a lattice with periodic boundary con- 
ditions. Since all stabilizers commute, the ground state 
degeneracy follows from counting the number of indepen- 
dent stabilizers! 35 * 36 * In a lattice of N squares (per flavor), 
there are 27V spins, N plaquettes, and N stars, so that 
the number of locally unconstrained spin— | degrees of 
freedom is 

N spin - N stab = (4) 

The non-trivial GSD follows from (spatially) global con- 
straints obeyed by the stabilizers. For each flavor, there 



are 2 global constraints: the product of all stars is the 
identity operator, as well as the product of all plaquettes, 

n^'=i (5) 

S 

l[Bf = l, (6) 

p 

because Pauli matrices act twice on each spin. The two 
unconstrained spin— I give 

\GSD ldeal2d \ = 2 2 -2 2 . (7) 

Importantly for our purpose, the GS manifold can be 
explicitly labeled using additional operators. There are 
two types (per flavor) of string operators^. Consider a 
path ("string") P/ along the /-flavored lattice edges, con- 
necting two star positions: all spins i on the links within 
the path are i £ Ff. Similarly, a dual path (string) Ff 
connects plaquette centers in the dual lattice of flavor /. 
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The string operators are then: 

4 d ( p /) = n 

'11"? 

ier f 



"2d 



(8) 
(9) 



If a string forms a closed loop, the string operator com- 
mutes with all stabilizers, and therefore with H, because 
a string without endpoints must share an even number of 
spins with any stabilizer. Further, a string shape P/P on 
the lattice can be deformed by multiplying the string op- 
erator T.2d /^2d by B/A stabilizers (see Fig. [I]). Such an 
operation preserves the dependence or independence of 
the string operator from the stabilizers. Then, Hidl'^2d 
is the trivial operator if P/P is contractible on the lat- 
tice in that way. The only non-contractible closed strings 
(loops) span the periodic system. These string loops are 
therefore independent from the set of stabilizers, and can 
be used to label the degenerate GS. A loop spanning 
the X (Y) direction has an odd number of intersections 
with a Si, string spanning Y (X), which causes their an- 
ticommutation. Due to their smooth deformations (using 
stabilizers), independent string loops are exhausted by 
topologically inequivalent ones. Any loop with multiple 
windings along the system can be deformed into a prod- 
uct of the elementary X and Y spanning ones. Therefore, 
a single flavor has two pairs of anticommuting I^^E^d 
(the X, Y and Y, X spanning pairs) , giving the expected 
total GSD from Eq. 0. 

We can now move on to a lattice with dislocations. 
Generally speaking, dislocations in 2d are pointlike, only 
of edge type, and are operationally obtained by remov- 
ing a finite line (the Volterra line) of lattice sites and 
repairing the lattice between the line's endpoints (see 
Fig. [I]). The endpoints mark two dislocations, transla- 
tional topological defected! We only consider disloca- 
tions with Burgers vector b = y/2, and therefore only 
remove single lines of sites lying along the x direction in 
Fig. [I] The latt ice is locally repaired seamlessly across 
the Volterra line away from its endpoints, and so the line 
itself is unphysicaP^, and only its endpoints matter. The 
Burgers vector topology says that any loop encircling a 
single dislocation will experience the "jump" due to the 
removed sites, even though it is impossible to say where 
the removed sites were. 

In the case of the 2dTC, our choice of b leads to the 
mixing of the two copies of the model (the two flavors 
/) explicitly through the local repair of stabilizers near 
the removed line (example of A s in Fig. [I]). However, 
strictly speaking, the global assignment of flavor becomes 
impossible due to the topological obstruction. We can 
however always assign flavors locally (starting from a de- 
sired point), and use it to build up arbitrary operators 
of fixed flavor step- by-step from the selected point, using 
the repaired local lattice links. The switch of flavors in- 
duced by a dislocation will always appear when we follow 
a string back to its starting point. Equivalently, and con- 



veniently, we instead retain the global flavor assignment 
of the original lattice together with an explicit choice of 
the removed site (Volterra) line. 

We can now repeat the counting of degrees of freedom 
from Eqs.Q, Around the dislocation positions the 
lattice is heavily distorted, and so the stabilizers there 
have to be carefully and explicitly reconstructed such 
that they commute with all others, as shown in Fig. [T] 
Along the Volterra line (shaded region of Fig[I]), each re- 
moved spin is also labeling a removed stabilizer. At a 
dislocation point however one stabilizer is lost (two orig- 
inal a,/3 stabilizers next to it merge into one). However, 
also two global stabilizer constraints are lost. Namely, 
the new stabilizers near the defects are still local but act 
on spins of both flavor, because of which now the flavor 
separated global constraints of Eq. ^ are merged 
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(10) 

(11) 



where the sets S, P, D s and D p exhaust all stabilizers 
in the dislocated lattice: original stabilizers and ones re- 
paired along the Volterra line are in the sets S and P 
(having well defined local flavor), while D s ,D p are the 
sets of heavily modified stabilizers near the dislocation 
cores. In a lattice with single dislocation pair therefore 
the changes in counting unconstrained degrees of freedom 
cancel, leaving the GSD intact. For k dislocation pairs, 
compared to ideal lattice, we get additional 2k degrees 
of freedom, but still only lose 2 global constraints. We 
therefore get additional 2k — 2 unconstrained degrees of 
freedom, predicting 



\GSD- 



2d 



\GSDide 



d2d\ 



-,2fe-2 



,k-l 



(12) 



for k > 1. Exact numerical lattice calculation as intro- 
duced in Sections |IV A[ |IVB| shows that 



\GSD 



2dnum 



2 4 • 4 k ~\ k > 1, 



(13) 



at least for k sC 3, and it also corroborates the inde- 
pendent string operators we will present in Fig. [2j The 
following discussion is based on a nice geometric inter- 
pretation in Ref 12. 

The neat feature of the 2dTC is that we can visualize 
the topology of the model manifold entirely. The 2dTC 
starts as two separate copies of the two-torus T 2 , one per 
flavor. The 71"! (T 2 ) = 1? gives the X and Y spanning 
non-contractible loops per flavored torus. 

Next we introduce a single dislocation pair. Because 
of the local flavor change upon crossing the Volterra line 
of removed sites (as explained above), we call this line 
the "a/3 line" . As an example, Fig. [T] shows an explicit 
string loop which winds around a dislocation. The string 
operator is constructed by definition, in a step-by-step 
manner: Starting from a star location and hopping to a 
neighboring star with which it shares a link, the spin site 



6 



on this link is acted on by a z and becomes a step in the 
path P of the string Eq. Q . The string is forced to wind 
(at least) twice around the defect, and in general an even 
number of times, because it must change flavor an even 
number of times to be able to close into a loop. 

Generally, knowing how strings change flavor upon 
crossing the a/3 line between dislocations, we see that 
the a/3 line represents a connection between the two fla- 
vored tori. This connection is a smooth tube, because 
we can show from the lattice model (e.g. using Fig. [I]) 
that a loop encircling the a/3 line (therefore both de- 
fects) can use the a/3 line to switch flavors, but cannot 
contract through it. Therefore the loop switches flavored 
tori by sliding along the a tube that connects them, see 
also Fig. ^p. 

Furthermore, the string loop winding twice around a 
single dislocation, shown as a string example in Fig. [T] is 
contractible (so giving the identity operator). This is a 
non-trivial fact which follows from the definition of heav- 
ily edited stabilizers near the dislocation site. Going to 
the continuum description, this loop can be contracted on 
the "connection tube" between the tori (this geometric 
interpretation was first introduced in Ref.38 as "worm- 
hole"), because it does not wind around the tube (and is 
therefore not a new type of non-contractible string loop) . 
Summarily, the k — 1 system can be topologically accu- 
rately described by a torus with two holes (the original 
holes of the two flavored tori), a manifold whose funda- 
mental group also has four generators and therefore the 
same GSD as the ideal system. 

Adding another pair of dislocations creates another 
tube connection between the two flavored tori, essentially 
adding another hole to the manifold, Fig. [2]d. This adds 
two anticommuting non-contractible string pairs, giving 
a GSD enhancement by factor 2 2 . One pair has T,2d en- 
circling the new dislocation pair and piercing both 
a/3 lines (therefore changing flavor twice); the other pair 
has T,2d and EE^d exchanged. Obviously the two strings 
of such pairs intersect once due to their flavor structures, 
see Fig. [2^,. 

For k ^ 1 dislocation pairs (Fig. [2ji) , the model lives 
on a manifold with 2 + (k — 1) holes, thus having genus 
g = k + 1. We name this two-manifold T 2 . Since 
iriiTg") — Z 2fl , there are 2g independent, say z^2d, strings. 
Another useful property of T 2 is that one can always 
choose the non-contractible paths in it such that they, in 
pairs, intersect exactly once. This ensures pairwise an- 
ticommutation with ^2d strings (which as a set do not 
topologically differ from the T, 2 d strings), and leads to a 
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2 2g = 2 2k+2 k ^ 1 



(14) 




b) 














FIG. 2: The GSD resolving non-contractible T,2d strings in 
the 2dTC (the H2d strings are topologically the same). The 
Volterra a/3 line is dashed blue, a/j3 flavored strings are 
straight/wavy black lines, dislocations are blue dots, a) Two 
dislocation pairs and the relevant strings, periodic boundary 
conditions are used, b) Another, equivalent, set of GSD re- 
solving strings, c) The underlying manifold of the model. 
Two tori are the two flavor copies of the system, and the 
dislocation induced a z strings from (a) are shown, d) The 
2 + (k — 1) genus manifold on which the k dislocation pair 
model lives, with the 2 new £2^ strings (of the type as in (b)) 
added per each new pair after the first. 



III. 



THE THREE-DIMENSIONAL DOUBLED 
TORIC CODE 



The 3dTC has been studied in detaiP^, but for com- 
pleteness we introduce all its necessary ingredients here 
(see Fig. [3]) as we describe its doubled variant. 

Spin— 2 degrees of freedom belonging to one flavor (say 
a) occupy the midpoints of edges in a cubic lattice, where 
the lattice sites s = (i,j,k) are labeled by integers, and 
lattice constant a = 1. Each lattice site s carries a sta- 
bilizer A s called a star operator (Fig. [3]) , which acts on 6 
spins neighboring s, i.e. on 6 spins positioned on lattice 
edges which share the site s: 



A* 



<J s+x/2 (J s 



-y/2 a s+z/2 s-x/2° s-y/2° 's-z/2' 



(15) 



where a l R is the i-th Pauli matrix acting on the spin at 
position R. Each square face of a cube is centered on 
p(ab) = s + (a + b) /2, for some s, and for a, b some pair 
from {£, y, £}; such a face lies in the ab plane and carries 
a plaquette operator Bp b acting on 4 neighboring spins of 

p(ab), i.e. acting on 4 spins on the lattice edges forming 
the square face (Fig. [3]): 



rjab,a 



p(ab)+d/2 a p(ab)+b/2 a p(ab)-a/2 a p(ab)-b/2' 



degeneracy, exactly in accord with Eqs.(13) and (12) 



In the following, it will mostly suffice to use the shorter 
label for the plaquette operator, where it is under- 
stood that the position p determines the plane ab in 
which the plaquette is lying. 

The f3 flavor of the model is created as a copy of the 
lattice translated by y/2 (Fig. [3]); /3 spin sites are then 
always at a position where there is an a stabilizer site, 
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FIG. 3: Three-dimensional doubled toric code (compare to 2d 
system in Fig. [Tj). Red/yellow dots are positions of spin— | 
of flavor a/P, and they differ by a translation by y/2. For 
convenience, we draw the cubic lattice having lattice spacing 
1/2 in all directions, so that spins, stars, plaquettes and local 
constraints, of both flavor, are all positioned on the sites of 
this "dense lattice". Flavor a//3 star operators (green/blue 
circles) are positioned on a//3 lattice sites and act through 
a x operators on their six neighboring spins (blue wavy lines 
from As). Flavor / = a, f3 plaquette operator Bp b '^ acts by 
o z on four spins of a / square lying in ab plane, centered on 
p. On the shown "dense lattice", a center p of a/j3 flavored 
square face is always positioned either: at an un-circled /3/a 
spin (example Bp V ' a shown), or at an empty position (center 
of j3/a cube, having no spin, nor star), example Bp X ' a . A 
closed string operator (example E°), forms the dashed loop 
P (example lying in YZ and ZX planes of a lattice) and 
acts by a z on spins in P (crossed- out). Multiplying E Q with 
shown Bp operators deforms the P loop. A center of a cube in 
the / flavored lattice is an empty position (no spin, nor star) 
in the shown "dense lattice", and it labels a local constraint 
Y\ Bp h,i = 1 for the 6 plaquettes p(ab) on the 6 sides of the 
given / flavored cube. 



and vice versa. Notice that every second XY layer of the 
3dTC model is reminiscent of a 2dTC lattice. 

The Hamiltonian of the 3dTC is simply given by the 
negative sum of all stabilizers in the lattice: 



H = — 




p(ab) 



(17) 



Setting if on a three-dimensional lattice with peri- 
odic boundary conditions, all stabilizers commute. The 



ground state degeneracy therefore fol lows from counting 
the number of independent stabilizers! 35 * 36 lEach center of 
cube in the / flavored lattice marks a local "cubic" con- 
straint, because the action of six / plaquette operators 
on the sides of that cube is the identity operator: 



ped 



1, 



(18) 



with Ci being one of the i = 1, . . . , N cubes in the / 
lattice, and p € C, being the six centers of squares on 
the surface of cube Cj. Per flavor, in a lattice of N cubes 
there are 3iV spins, 3iV plaquettes, and N stars, so that 
the number of locally unconstrained spin— | degrees of 
freedom is 



N, 



spin 



^indep.loc. constraints) — 



3N - (47V - N) = 0, (19) 



balancing as in the 2d case of previous section. However, 
only N — 1 local cubic constraints from Eq. ( 18 1 are in- 

-9 

1 



dependent, since due to 



1 the product of N 
local constraints gives exactly the iV-th constraint in the 
system with periodic boundary conditions. This removes 
one stabilizer constraint. On the other hand, the product 
of all stars in the lattice gives identity: 



(20) 



adding one global constraint (per flavor). The counting 
balance of Eq. (19) therefore still seems to remain valid. 

However, there are additional global constraints 
obeyed by the stabilizers, from which the non-trivial GSD 
follows. The stars are "volume-like" , since multiplying 
many adjacent stars leaves only a surface of spins acted 
on by a x . The plaquettes are "surface-like", since mul- 
tiplying many adjacent plaquettes leaves only a closed 
loop of lattice edges on the boundary, on which a z oper- 
ators act (e.g. the loop operator in Fig. [3] is equal to the 
product of several adjacent plaquette operators). We al- 
ready counted the global star constraint obtained by fill- 
ing the system volume by multiplying all stars (thereby 
shrinking its surface acted on by a x operators to noth- 
ing). The periodic system has a shape of the three-torus, 
and we can multiply plaquettes so as to create three sur- 
faces spanning the system's XY, YZ or ZX planes, such 
that they have no edges and so no a z action. This gives 
three global constraints. Each surface can be arbitrarily 
deformed by putting additional plaquettes into the prod- 
uct. But, demanding that the surface stays boundary- less 
(and therefore equal to 1 operator), the elementary defor- 
mation must contain the product of plaquettes forming 
the surface of an entire cube. The product of plaque- 
ttes on a cube surface is locally constrained (Eq. (18)), 
so the global constraint given by the deformed surface 
is, up to local constraints, equivalent to the global con- 
straint given by the undeformed surface. The three non- 
contractible (and boundary-less) surfaces spanning the 



three-torus therefore give the only inequivalent global 
constraints which are independent from the local con- 
straints. The existence of three spanning surfaces follows 
from 7Ti(T 3 ) = T? . Adding up two flavors, one gets 



\GSD id , 



ealSd 



= 2 b 



(21) 



on the three-torus, i.e. on a lattice with periodic bound- 
ary conditions. 

The operators which resolve the ground state manifold 
are represented by non-contractible str ing lo ops (Fig. [3]) 
and non-contractible closed membrane d 26 l 34 l A string 



E'flv) - n < 



(22) 



is a product of operators on spins positioned along a path 
P/ that follows the edges of / € {a, f3} flavored cubic 
lattice, starting and ending on a lattice site s, on which 
there is a star A{ (Fig. kjl). The string can obviously 
be deformed by multiplication with plaquettes, and is a 
straightforward generalization of 2d strings X^d- If the 
string is closed, it commutes with all the stars too, as 
every star acting on an edge belonging to the string will 
also act on another of string's edges due to the string's 
continuity. A membrane 



r'(s/) = n < 



(23) 



is defined on the dual lattice, and is built from opera- 
tors acting on a single spin per dual cube face (see il- 
lustration in Ref.34). The membranes can be deformed 
when multiplied by stars, which are in the centers of dual 
cubes. (This deformation is geometrically the same as 
when the global surface constraints on plaquettes are de- 
formed by multiplication with local cubic constraints.) 
Since the membranes are closed (boundary-less), they 
commute with all plaquettes, and are generalizations of 
the 2d strings hopping on the dual square lattice. 

By construction, a string shares exactly one spin with 
a membrane at their intersection point, which leads to 
their anticommutation in case the string pierces the mem- 
brane odd number of times. All closed membranes must 
be pierced even number of times by any string loop, ex- 
cept in the case of non-contractible objects which exploit 
the three-torus topology. On the dual cubic lattice, one 
can construct only the three independent closed mem- 
branes Tf(XY),Tf{YZ),Tl(ZX). Each is pierced ex- 
actly once by £^(2T), E^(A), Y,f (Y) (strings spanning the 
three-torus), respectively, giving 3 anticommuting pairs 
per flavor. These strings and membranes commute with 
stabilizers and the Hamiltonian, so each anticommuting 
pair gives a two- valued quantum number labeling the de- 
generate ground states. Due to flavor / there are in total 
6 anticommuting pairs, giving the \GSDid e ai3d\ = 2 6 . 




FIG. 4: Edge dislocation of an ZX loop (compare to Fig. [TJ. 
Black dots are removed sites, forming the a/3 surface shaded 
grey. The dislocation line is represented by the outmost array 
of black sites, marked by the blue line. The shown star opera- 
tor A s is repaired by just gluing the neighboring sites together 
in the local neighbor network, and that is how the flavor jump 
occurs as a string operator which pierces the a/3 plane is con- 
structed step-by-step (compare to Fig. [I]). The exception to 
the seamless repair are stabilizers that contain spins on the 
dislocation lines, for which the explicit repair "surgery", such 
that all stabilizers commute, is needed, as shown in Fig. |16| 



A. Dislocation loops 

In this subsection we present the details of creating dis- 
locations in the 3dTC model, which is a generalization of 
the 2d case from previous section. Due to our choice of 
displacement of flavored lattices along y, we only consider 
the nontrivial dislocations with b = y/2. The dislocation 
loops in ZX plane are therefore purely edge (we will call 
them "edge dislocation loops"), while the XY and YZ 
dislocation loops are of screw type along segments par- 
allel to y segments, and edge otherwise (these we simply 
call "screw dislocation loops" ) . 
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To create an edge dislocation loop^ 1 , we remove a sur- 
face of sites in a ZX plane, as in example of Fig. [5] (see 
Figs. 16|17 for details). Locally, the lattice, as well as 
the stabilizers, depend on the nearest neighbor network, 
which can be repaired by connecting two sites on either 
side, along y direction, of a removed site. This reparation 
is locally indistinguishable from the original lattice, but 
fails at the dislocation line formed by the outermost ar- 
ray of removed sites, as dictated by the defect topology. 
The non-trivial reparation of stabilizers around the dis- 
location line is presented explicitly in Appendix [El Once 
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"repaired" , all stabilizers in the lattice commute, and we 



use their negative sum for the Hamiltonian as in Eq. ( 17 1 



tion for the dislocation influence on global constraints 
(see Section IV C I. 



The removal of a single layer of sites introduces a 
"skip" in the usual alternation of flavors along y, thereby 
inducing the flavor mixing. To be precise, the plane sur- 
face of removed sites is the Volterra surface, which is un- 
physical, i.e. could be arbitrarily placed in the lattice as 
long as it is bounded by the fixed dislocation line^. This 
also means that the global assignment of flavor is impos- 
sible in presence of dislocations. However, we can always 
locally assign flavors in the lattice. Most importantly, 
following a string that starts out locally as flavor a, fol- 
lowing the local neighbor network, if it winds around a 
dislocation line and comes back near the start position, 
the "skip" will be felt (although impossible to say where 
it happened) and the string will not close on itself as it 
will locally be of j3 flavor. (We need to perform one more 
wind around a dislocation line to be able to close the 
string.) 

We can alternatively keep the original global flavor as- 
signment, but then also keep the surface of removed sites 
as a fixed Volterra choice, for the simplest way for book- 
keeping. Because of this, we will also call the chosen 
Volterra surface "the a/3 surface" as the flavor change 
will occur upon crossing it. This procedure only makes 
sense when considering closed strings and membranes 
which conserve the parity of number of times they cross 
the a/3 surface that spans a dislocation loop. Since our 
goal is analysis of GSD, the open strings and membranes 
carrying excitations on their endpoints/boundaries^will 
never be involved. 

For the creation of screw dislocation loops, we need to 
perform ay/2 translation which is parallel to the chosen 
"screw surface" tha t spans the dislocation loop (no sites 
are removed), Figs. 18|19 This translation in the local 
neighbor network also introduces the "skip" between fla- 
vors (along y), just as for the edge loops. All remarks 
concerning the a/3 surface remain as in the edge disloca- 
tion case. 

The preservation of the local neighbor network pre- 
serves the commutation of stabilizers involved in the dis- 
location surface. However, the stabilizers involving spins 
on the dislocation line need to be repaired, removed or 
replaced in a non-trivial way. We explicitly construct the 
new stars and plaquettes for all types of dislocation lines 
and their corners. We leave the presentation of details 
for the Appendix, Figs. [I6p0l The main outcome is that 
the recipe for lattice "surgery" produces the new local 
stars and plaquettes as well as local cubic constraints in 
such a way that the local counting of unconstrained de- 



grees of freedom remains just as in Eq. ( 19 ). This is true 



also for linked dislocation loops (when two a/3 surfaces 
intersect). However, the dislocations change the global 
properties of the model. In our exact numerical calcu- 
lations (Section IV A) we corroborate the local balance 
of Eq. ( 19 ) by showing that the GSD in particular does 



not depend on the sizes and shapes of simple dislocation 
loops in the lattice. We further corroborate the predic- 



IV. RESULTS 

This section has three parts: We start by present- 
ing our exact numerical calculation, which includes the 
GSD for various dislocation configurations, as well as 
the demonstration of independence of relevant string and 
membrane operators; next we present the analytical ar- 
guments that lead us to the particular choice of the 
strings and membranes to analyze, and we also give the 
prediction of the GSD; and finally we close with all re- 
sults for the much simpler two-dimensional model. 



A. Exact numerical results 

We analyze 3dTC on lattices with various dislocation 
configurations, as introduced above and in detail in the 
Appendix (especially captions of Figs. I6p0 |. The lat- 
tices have variable sizes, typically about 250, and up to 
720 cubes, which are sufficient to study all presented sce- 
narios. 

We find the GSD directly using a quantum informa- 
tion theorerri- 35 ! 36 ^ Counting the number of spins in the 
(dislocated) lattice, and subtracting the number of in- 
dependent stabilizers in that lattice gives the number of 
"unconstrained" spins; the GSD then equals the size of 
the spin (which is 2 for spin— 1/2) raised to the number 
of unconstrained spins. 

The procedure for finding the number of independent 
stabilizers (stars and plaquettes separately, but including 
both flavors) we used is also direct. For example, con- 
sider first the stars, with N star being their total number. 
Each stabilizer is represented by the list of spins on which 
it acts. Since this star operator is just o~ x or 1 on any 
particular spin site, we just need to assign 1 (0) if the 
o~ x acts odd (even) number of times on the given spin. 
Therefore, on each spin, the operator is represented by 
the additive Z 2 algebra. There is a total of N s spins in 
the lattice, so we next form an N star xN s matrix, the i-th 
row of length N s filled with 1 for the spins on which the 
i-th star acts, and filled with on all other spins. Each 
star contributes a row to the matrix, as a candidate for an 
independent row under the Z 2 algebra. Finally, the ma- 
trix is converted to the row reduced echelon form under 
the algebra, and its rank found (or number of zero rows 
counted), which is the number of independent operators. 
The same procedure is used for the plaquettes. 

We numerically studied up to four dislocation loops of 
varying sizes and shapes each. The linking did not influ- 
ence the outcome, which only depends on the number of 
loops; among up to four loops, we could check: no link- 
ing; one linked pair; two separated linked pairs; three 
linked loops; two parallel screw loops, each linked with 
the two parallel edge loops. 
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We get \GSD i( ieai3d\ — 2 6 in the ideal lattice, as ex- 
pected from Eq. (21), and \GSD\ = 2 6 • 2 k ~ l for k dislo- 



cation loops in the cases k 6 {1, 2, 3, 4}. 

A special case occurs when a dislocation line closes 
into a loop by spanning the periodic system (these lines 
need to occur in pairs). For one pair we get \GSD\ = 2 s , 
less than in the ideal lattice. The GSD does not change 
further by introducing spanning dislocation pairs in an 
orthogonal direction. However, we checked that upon in- 
troducing k £ {2, 3, 4} parallel spanning pairs, the count- 
ing follows a trend similar to the 2dTC case j Cff-Pl = 



for 



2 • 2 2k+2 = 2\GSD 2d \ (see ends of Sections |iVD||iVC 
details). 

We performed the same numerical analysis on the 
2dTC model, where we considered up to three disloca- 
tion pairs. 

Finally, we constructed a 3dTC lattice with a dislo- 
cation line knotted into the Trefoil knot, as shown in 
Fig. [8^,. As our analytical arguments in the next section 
show, this knotted dislocation behaves as a simple dislo- 
cation loop, and indeed in the exact numerical calculation 
it gives the corresponding | GS'Z) | = 2 6 . 

The summary is presented in Table [I] 



Configuration 


GSD 


3dTC 




Ideal 


2 6 


k ^ 1 loops (independent of linking, knotting) 


2 6 . 2 k-l 


1 loop 


2 6 


1 spanning loop 


2 5 


2 orthogonal spanning loops 


2 5 


k ^ 2 parallel spanning loops 




Trefoil knot dislocation 


2 6 


2dTC 




Ideal 


2 4 


k ^ 1 dislocation pairs 


2 4 . 4 fc-l 



TABLE I: GSD for k < 4 dislocation loops in the 3dTC, and 
k ^ 3 dislocation pairs in the 2dTC, found in exact numerical 
calculations. 



B. Exact numerical analysis of stabilizers 

Using the above method, we could also check directly 
the independence of various strings X* , membranes E* , 
and global plaquette constraints. Namely the strings can 
be deformed by multiplication with plaquette operators, 
and therefore two equivalent strings will not appear as 
independent additions to the set of all plaquettes. The 
same applies for the set of E/ and A? . In particular, 
any contractible X (S) will also not appear independent 
when added to the set of B (A). 



k = 4 loops, that all the predicted anticommuting non- 
contractible strings and membranes are indeed indepen- 
dent within the set of plaquettes and stars, respectively. 

We delay the detailed presentation of practical con- 
struction of membranes on the lattice with dislocations 
to Appendix [5] because the analytical proof that follows 
is presented through the string perspective. 



C. Continuum strings 

In this subsection we analyze the behavior of GSD by 
counting non-contractiblc strings, and give the prediction 
for GSD of unknotted dislocation loops with allowed link- 
ing. 

Our present goal is to find the set of independent 
non-contractible string operators. Such strings neces- 
sarily pairwise anticommute with a set of independent 
non-contractible membranes (string piercing membrane 
odd number of times), thereby enumerating the GSD. To 
streamline the presentation, we therefore only count the 
strings, and delay the (in itself illuminating) analysis of 
membranes to Appendix |B| 

To this end, we have derived, by direct analysis of all 
situations on the 3dTC lattice (Figs. 16p0 |, a set of gen- 
eral rules that the strings obey as they are deformed and 
moved through the dislocated lattice by the action of pla- 
quettes. We of course use the plaquettes obtained by the 
proper "surgery" in the dislocated lattice (Append ix |E|. 
Having in mind the detailed discussion in Section [HlAl 
we fix the Volterra a/3 surface of each dislocation loop, 
and retain the global flavor assignment from the original 
ideal lattice. 

The fundamental set of lattice rules we derived for the 
behavior of strings is shown in Fig. [5] Intuitively, the 
string is an elastic reconnectable "rubber band" which 
switches flavor only upon crossing through an a/3 surface. 
The string cannot pass through a dislocation line, but can 
be contracted if it winds an even number (w £ 2!L) of 
times around the same dislocation line. The last property 
is crucial, but not intuitive. Deriving it demands the 
precise use of the stabilizers obtained by "surgery" on 
the lattice at the dislocation line. (A similar property 
holds in the 2dTC case, see Fig.[l}) 

We analyze the consequences of these derived rules on 
two simple but insightful examples, obtaining outcomes 
that match exact numerical lattice results: 



Using the analytical predictions (Section IV C), we 



could indeed check numerically, in lattices with up to 



• A single dislocation loop does not change the GSD. 
First, we note that the standard strings which span 
the system in 3d are still independent, and free 
to roam through space avoiding the finite dislo- 
cation loop. However, every new candidate non- 
contractible string loop must pierce the a/3 sur- 
face at least twice (or generally an even number of 
times): the string has to change flavor even num- 
ber of times along its length to be able to close 
into a loop. But then it also winds an even num- 
ber of times around the dislocation loop line, so 
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ing operator pair when an appropriate membrane is 
considered, leading to GSD enhancement by factor 
of 2). The new candidate non-contractible string 
can now exist because there are two separate a/3 
surfaces. The candidate pierces both of them once, 
therefore winding only once around each disloca- 
tion loop line. This string is not contractible, as we 
cannot change the parity of times it pierces any a/3 
surface, neither by deformation, nor by attempted 
(but forbidden) passage of the string across a dis- 
location line. 

In the two dislocation loop case, it is amusing to con- 
sider a loop that does not wind around either dislocation 
line, yet seems "stuck" because it is linked with them. 
We consider such a Borromean ring configuration of two 
dislocation loops and one string loop, Fig. [SJs. Because 
the string pierces through each a/3 surface an even num- 
ber of times, the string is contractible. We also check 
this numerically on the lattice by showing that such a 
string is not independent from the set of plaquettes on 
this lattice. 



FIG. 5: Elementary allowed operations on string operators in 
presence of dislocation loops, derived from their properties on 
the lattice. Straight and wavy black lines represent strings of 
two flavors, / = a, /3, blue lines are dislocation lines, and the 



dotted surface is the Volterra a/3 surface (see Section III A I. 
Dots mark the piercing of string through a/3 surface, (a) Re- 
connection of two string segments, (b) Local string flavor 
change using dislocation; a top view detailing the process is 
also shown, (c) Local passage of string through dislocation 
line, (d) String winding w — 2 times around dislocation line 
is contractible; only even w is possible due to flavor continu- 
ity, all contractible. (e) An example of contracting a string, 
using rules (a) and (d), for a string (black line, but its fla- 
vor marked only after second subfigure) which is non-trivially 
linked with two dislocation loops (blue and red) in a Bor- 
romean ring configuration. Green symbols mark the points of 
reconnection (rule (a)). The string is therefore not indepen- 
dent from plaquette operators, i.e. can be written as some 
product of those. 



it is contractible according to the special derived 
rule in Fig. [5ji. Imagine an alternative candidate 
string that spans the system twice in one direction, 
piercing the a/3 surface twice along the way. We 
can combine it with both flavored (and purportedly 
independent) standard system-spanning strings in 
the same direction, to reconnect them into non- 
system-spanning loops, of which one of either fla- 
vor will be far from the dislocation and trivially 
contractible; the leftover candidate is just the type 
we considered (winding locally around the disloca- 
tion). Other combined possibilities can also clearly 
be reduced to these cases. 

Two separate dislocation loops add one non- 
contractible string (and therefore one anticommut- 



D. Proof of GSD and non-contractible strings for k 
unknotted dislocation loops 

One can now generalize to k separate dislocation loops. 
Here we consider simple dislocation loops, having ori- 
ented, non-self-intersecting a/3 surfaces. It is simple to 
see that the addition of a new separate dislocation to the 
system adds one independent string, increasing the GSD 
by factor 2. Focus on a candidate for a new independent 
string; it must pierce the new a/3 surface an odd num- 
ber of times, since we will now show that otherwise it 
could be deformed until completely separated from the 
new dislocation loop, becoming therefore equivalent to 
some string in the old set of strings. The odd num- 
ber In + 1 (n ^ 0) of piercings of the new a/3 surface 
is consequently equivalent to a single piercing. To show 
this, consider a candidate string with an even number 2n, 
n ^ 1, of piercing points through the new dislocation's 
a/3 surface. Each of 2n piercing points has two string 
segments of opposite flavor emanating to the two sides of 
a/3 surface. We can deform the string in the vicinity of 
the new dislocation loop to reconnect the even number 
2n of emanating segments of cither color until we obtain 
a collection of n disconnected loops, each containing ex- 
actly two piercing points. These loops with two piercing 
points are contractible either trivially, or by the nontriv- 
ial rule if they wind w = 2 times around the dislocation 
line (see also Fig. [8| . All the 2n piercing points are thus 
removed, and the string remains completely separated 
from the new dislocation loop. 

So the set of candidates for new non-contractible 
strings boils down to strings which pierce the new dis- 
location loop once. These are however all equivalent. 
Imagine there are two which are independent. They both 
pierce the new loop once, and we can deform them to 
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coincide in the region around the piercing point. By re- 
connecting them at the points where they deviate from 
each other, they just cancel on the piece containing the 
piercing point (where they coincided), while the remain- 
ing piece is a string that is completely separate from the 
new dislocation loop. It therefore must be dependent 
(a combination of) the old strings. The two new can- 
didates differ only by an old string, and are therefore 
not independent of each other in the enlarged set of all 
non-contractible strings. 

The linking of dislocation line loops surprisingly leaves 
the situation unchanged. Consider two linked dislocation 
loops: their a/3 surfaces intersect along some finite line 
(see Fig. 14 3 for illustration). However, a string pierc- 
ing each of the intersecting a/3 surfaces once cannot be 
contracted. Namely, the line of intersection of the a/3 
surfaces does not allow the passing of the piercing points 
from one surface to the other. Therefore, the linked loops 
behave as separated loops, at least in terms of the topol- 
ogy of strings, membranes, and the GSD. 

Adding the extra factor 2 to the GSD for each dislo- 
cation loop after the second one, gives the numerically 
corroborated 



\GSD\ 



2 6 • 2 k 1 for k ^ 1 (possibly linked) loops. 

Finally, a dislocation line pair which spans the system 
along some direction (say Z) warrants a note. The GSD 
is lowered — there are no new independent string candi- 
dates, but instead one non-contractible string spanning 
the system along Z is lost. Namely, a string spanning 
Z direction can switch flavors by passing in its entirety 
through the a/3 surface which also spans the Z direc- 
tion. Therefore the two flavors for the Z spanning string 
are not independent anymore. The system-spanning XY 
membranes of two flavors also merge into one when such 
dislocations are present. For an array of such disloca- 
tions, one can make a mapping to the 2dTC, as also 
discussed from the membrane viewpoint in Appendix [B| 



E. Dislocation knotting 

This subsection analyzes the case of knotted disloca- 
tion loops. Since we dealt with linking previously, it is 
enough to study a single knotted dislocation line loop. 

Before we show why dislocation knotting in general 
does not influence the GSD value, we need to focus on 
the nature of the Volterra a/3 surface, which becomes 
non-trivial for a dislocation line twisted into a knot. 

There are several ways of setting a general Volterra 
a/3 surface which spans a knotted loop. By the phys- 
ical nature of the Burgers vector, the a/3 surface must 
be orientable, so we disregard the nonorientable option 
(for the Trefoil knot which is not carrying a Burgers vec- 
tor c harg e, the triple twisted Mobius strip could be an 
optiorP^l, see Fig. Up). However, the surface can still be 
self-intersecting, and we discuss such cases at the end of 
this subsection. 



We first consider the non-self-intersecting option for 
the a/3 surface. Remarkably, such an orientable, con- 
nected, non-self-intersecting surface span ning an arbi- 
trary knot can always be constructe d 39 * 40 *, and is called 
the Seifert surface. The explicit construction is based on 
a 2d projection of the knot, which is cut into discs, that 
are finally connected by twisted strips (see Fig. [7^, for 
the Trefoil knot example). The resulting surface is non- 
unique and even has a non-unique genus (bounded from 
below), but it still satisfies the important criteria from 
above. Its (single) edge is topologically equivalent to the 
knot one started from. One can in principle construct 
such an a/3 surface even on the lattice, because it is pos- 
sible to gradually twist the orientation (say in the XY 
plane) of a strip of selected lattice sites, as one follows 
the strip along its length (say Z axis). 

Certainly, we can use the Seifert a/3 surface to prove 
the absence of non-contractible string candidates. The 
proof is straightforward, as all the elements have already 



been introduced in Section IV D The key properties are 
that: 

PI There is a single, orientable and connected a/3 sur- 
face. 

P2 There is a single dislocation line (the edge of the a/3 
surface). 

The main steps of this proof are presented in Fig. [6] for 
the example of a dislocaiton making a Trefoil knot. 

We start from a candidate string that pierces the sur- 
face, necessarily even number of times, say m — 2n times, 
due to string continuity and flavor changes. Using prop- 
erty PI, one can deform the string to bring all 2n piercing 
points next to each other on the a/3 surface, and close 
to some chosen point R on the dislocation line, by us- 
ing only smooth string deformations. (The deformation 
might make the string even more tangled-up, but this 
does not matter.) One then locally (near the point R on 
dislocation line) reconnects the bundle of string segments 
emanating from the piercing points (2n points, giving 2n 
segments per flavor). The reconnection, in the sense of 
Fig. [5^,, breaks up the string loop into a set of discon- 
nected string loops. Only three kinds of strings loops 
can appear: 1) Small string loops that wind around the 
dislocation line near the point R; 2) Small string loops 
that pierce the a/3 surface near point R, but do not wind 
around the dislocation line; and 3) Potentially long and 
tangled-up string loops that have no contact with the 
a/3 surface. The strings of type (1) and (2) are immedi- 
ately contractible, the former using the special rule from 
Fig.[5]i, and the latter trivially. The strings of type (3) we 
call "w = string loops" , since they do not pierce the a/3 
surface, and therefore do not wind around a dislocation 
line in the sense of Fig. [5ji. 

Our goal is to show that an arbitrary w — string 
loop is also contractible, which we do using a procedure 
demonstrated in Fig. [6j Focusing on one such string, it 
can either be completely detached from the knot, and 
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FIG. 6: Contracting a string that is tangled with a dislocation which forms a Trefoil knot, where the string does not pierce 
the a/3 surface. This is an example of contracting a w = string loop, needed for the proof in Section IV E a) 1. A w — 
string, to be reconnected at the red symbols. 2. It reconnects into four elementary string loops, cycles which pass through 
neighboring holes in the a/3 surface (and do not pierce it). 3. Focusing on one (marked in previous step) of those cycles, one 
uses the single dislocation line to reconnect it into two non-piercing string loops. 4. A string loop changes flavor when pushed 
through the surface. 5. The doubled loops can again be reconnected into elementary string loops (in this example no need, as 
it is already such). 6. Every elementary loop is doubled, and these can be contracted locally using the local a/3 surface, b) 
Zoom-in on typical operations needed for steps 3,4 and 5. 



therefore trivially contractible, or it can be entangled 
with the dislocation line knot by winding through the 
holes in the Seifert surface. The string is of single flavor 
because it does not pierce the af3 surface. This makes it 
simple to reconnect the string loop until we obtain a set of 
"elementary" string loops, namely the string loops that 
are elementary cycles which enter and exit the Seifert sur- 
face through a pair of neighboring holes in the surface. 
This is always possible to do according to the properties 
of Seifert surface in PI. An example of an elementary 
loop would be a string winding around one of the three 
twisted strips in the canonical Seifert surface of the Tre- 



foil knot, and such a string example is marked in Fig. [6^2. 
Fig. [6^,1,2 show the entire operation. 

Let us therefore focus on a single elementary string 
loop, and show that it is contractible. The elementary 
loop is tightened close to the Seifert aft surface, and by 
its definition, it can be seen to wind around two short 
segments of dislocation line which span a short piece of 
Seifert surface between them. These two dislocation line 
segments are centered on points A and B of the disloca- 
tion line, respectively. Since we will use these objects in 
the vicinity of points A and B, the two short segments of 
the dislocation line can be considered straight and paral- 
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in Fig. [8^, we made the opposite choice. The outcome 
is that the Trefoil knot constructed on the lattice using 
a self-intersecting Volterra surface does not change the 

GSD: \GSDTrefoil\ = \GSDideal3d\ = \G SDk=\disl. loops] ■ 

It seems therefore that the self-intersecting case be- 
haves just as we expect from the non-self-intersecting 



FIG. 7: The Volterra ("a/3") surface of a dislocation line 
twisted into a Trefoil knot, a) A non-self-intersecting and 
orientable a/3 surface always exists, e.g. the Seifert surface 
(left and center, two colors for two sides of the surface), b) 
A non-orientable surface spanning a knot can also sometimes 
be found, and here it is the triple twisted Mobius strip. See 
Fig. [8^ for an implementation of the Trefoil knot through a 
self-intersecting a/3 surface. 



lei, and the piece of Seifert surface they span locally flat, 
as shown in Fig. [6|d1|6£i2. 

We next manipulate two halves of the elementary 
string loop as shown in Fig. [6^,3j6]3l-4: one half stays 
hovering around the surface at point A; the other half 
of elementary string loop we drag all the way along the 
dislocation line starting from point B until we reach the 
point A, which is always possible due to property P2. 
What is created by the dragging is two parallel string 
segments connecting the half at A to the dragged half 
(now also at A). The locally parallel string segments are 
positioned tightly below and above the surface (which is 
orientable!). This is illustrated in Fig. [6^,3, b2,b3. These 
two segments we call d and u, respectively, and by con- 
struction there must be a piece of Seifert surface between 
them all the way along their length. At point A we can 
now reconnect the string loop, so that the "stay-at-A 
half" and the "drag-from-S half" cancel, while each u 
and d connect into a complete loop, see Fig.[6]33,4. The u 
and d are geometrically almost identical: they are offset 
by a short distance having the Seifert surface between 
them. Since the surface is orientable, we can push d 
through it, and according to standard rules end up with 
two copies of the u string loop, carrying opposite flavors. 
These twin strings are again just w = string loops and 
we can decompose them into elementary cycles, just as 
we did with the original w = string loop. The ele- 
mentary string loops that we obtain now consist of two 
copies each, having different flavors. The earlier situa- 
tion in Fig. now becomes the one in Fig. [6]d6. The 
elementary loops can make use of the flavor flip property 
(Fig. [5]d), and therefore we end up with two copies of 
the string having the same flavor, which just cancel each 
other. 

This completes the argument that any string loop con- 
structed in presence of a knotted dislocation line loop is 
contractible, and therefore can be written as some prod- 
uct of plaquette operators. 

So far in this subsection we focused on non-self- 
intersecting a/3 surfaces. However, for the efficient imple- 
mentation in the lattice code (Section IV A I , as sketched 



case. Recall from discussion in Section IV C and Fig. 14 
that from the lattice model we derived the rule saying 
that a string is not allowed to move freely across an in- 
tersection line of two pieces of a/3 surface. This situation 
seems to invalidate the first step of our proof, where we 
demanded free movement of the piercing point of string 
on the a/3 surface. In particular, this raises the ques- 
tion about the seemingly non-trivial candidate for a non- 
contractible string, which we construct such that: 1) It 
is entangled with the dislocation line; 2) Does not pierce 
the a/3 surface; and 3) Seems "stuck" (non-contractible) 
because the intersection line of two pieces of a/3 surface 
is in its the way. An example relevant for the Trefoil knot 
on the lattice we used is shown in Fig. [8] However, such a 
string loop can be easily locally contracted irrespective of 
the a/3 self-intersection line, as demonstrated in Fig. [SJd 
using the above general recipes. 



V. EMERGENT NON-ABELIAN GAUGE 
THEORY IN DISLOCATION MELTED DOUBLED 
TORIC CODE IN 2D 



In this section we consider the situation in which dislo- 
cations of the doubled Z„ toric code become dynamical, 
and show that the resulting phase of matter is governed 
by a non-Abelian gauge theory. (The Z„ toric code lattice 
model and its GSD, without melting, are described in Ap- 
pendix[Aj) The physical mechanism behind the transition 
is melting of the lattice through proliferation of "double 
dislocations" (lattice dislocations of length |b| — a). Con- 
densation of point-like electric and/or magnetic charges 
can lead to states with different topological order^H, how- 
ever in the present context the condensed double dislo- 
cation (point-like in 2d and loop-like in 3d) itself is not 
carrying any charge or flux in the theory. 

Our first step is to describe how the dislocations |6| = 
a/2 obtain dynamics by lattice melting, conjecturing that 
the resulting state is described by a G n = (Z„ x Z n ) xi Z 2 
gauge theory. We give a physical argument for the melted 
state GSD on a surface with arbitrary genus g, consider- 
ing first the case of torus. We next obtain the exact same 
GSD values for states described by gauge group G n on 
genus g surfaces. This implies that the quasiparticle con- 
tents also match, giving strong evidence that the melted 
state is really described by G n topological order. As a 
special case, this implies that the melted phase of the 
doubled Z 2 toric code is described by the non-Abelian 
group D4. 
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a) 




FIG. 8: Self-intersecting a/3 surface and string contraction, a) 
Self-intersection (along green line) of the a/3 surface (dashed 
area) in the Trefoil knot configuration we used on the lattice in 
Section [IV A| Some parts of the surface are shaded for better 
visibility, b) Contracting a string tangled with an a/3 surface 
that has a self-intersection (along green line). This case arises 
in the knot implementation from (a). The two shown surface 
pieces (shaded green and violet) could also be disconnected, 
since the string contraction is done locally. Piercing points 
of the string are marked with dots of corresponding color. 
The procedure hinges on moving piercing points close to each 
other and close to dislocation lines, but without passing the 
piercing point from one surface to the other, which is a for- 
bidden deformation. The string is then reconnected until it 
gives disconnected pieces that wind w — 2 times around a 
single dislocation line. 



A. Dynamic dislocations and physical argument for 
torus GSD 

Before considering the melted phase, we need to de- 
scribe the "twist" effect introduced by dislocations. Con- 
sider first a periodic 2d lattice, which is topologically 
equivalent to a torus. Before any melting, a single 
dislocation — anti-dislocation pair can be created, moved 
around the periodic system, and annihilated, as shown 
in Fig. |9^,b, specifically using the glide motion along 
b ~ y. This operation leaves a "Z2 twist loop" (blue 
dashed line in Fig. ^p,c) spanning the y direction of the 
system, along which all x rows are shifted by a/2 along 
y, and across which a flavor change is imposed. This im- 
plies that strings winding around x direction are forced 
to change flavor twice before being able to close, while 
strings winding along y change flavor when crossing the 
twist (Fig. [9b). The two flavors have merged, and the 
system is described by a single copy of the toric code on 




FIG. 9: Dynamic dislocations introduce Z2 twist in the dou- 
bled toric code. In unmelted lattice (red and yellow dots 
schematically depict the two vertically displaced flavors of 
spins in the doubled model) dislocations can (a) nucleate and 
(b) glide until annihilation, to create a "twist loop" (dashed 
blue line) across which all quasiparticles (strings) switch fla- 
vor, (c) In the continuum description of periodic system, 
string loops (i.e. quasiparticle tuneling operators) of both 
flavor (solid and wavy line) are shown. The "twist" replaces 
two copies of a torus by a single torus copy with L x — >■ 2L X . 
The two strings along y are actually equivalent due to the 
twist, (d) When the lattice melts, distinct topological ground 
state sectors are labeled by the twists they contain. In a genus 
g surface all sectors containing any twists topologically map 
to the shown case. Since dislocations delocalize in the melted 
phase, physical ground states labeled by strings tied to any 
of g — 1 holes have to be symmetrized in flavor. The physical 
states labeled by strings tied to the g-th hole are automati- 
cally symmetrized in flavor since the twist replaces the two 
flavor copies by a single one, as in (c). 



a torus of size 2L X x Ly. 

Next we consider the effect of melting the lattice 
through double-dislocations. Before melting, the fla- 
vor exchange symmetry is realized through ay/2 lattice 
translations. The double-dislocation = a) melted 
phase restores full translational symmetry, so that also 
flavors become locally indistinguishable. Since we are 
not interested in liquid crystal phases, a lattice with dis- 
locations is analogous to a U(l) superfluid with vortices. 
Proliferation of double-vortices leads to an insulator de- 
scribed by Z2 gauge theory, where the original single vor- 
tex becomes the deconfined gauge charge (and vortices in 
the double- vortex condensate are deconfined fluxes). By 
analogy, the original |6| = a/2 dislocations are deconfined 
gauge excitations in the melted phase. In such a state, 
we can say that the flavor-changing Volterra lines con- 
necting \b\ = a/2 dislocation — anti-dislocation pairs also 
permeate the system, again concluding that it is unphys- 
ical to locally distinguish the two flavors. Physical states 
are only the ones symmetric in flavor, i.e. the flavor Z2 
symmetry is gauged. 

The physical states can still be distinguished according 
to the Z2 twist loops they contain. Namely, a state with 
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a twist loop created by a non-contractible tunneling of a 
dislocation — anti-dislocation pair is topologically distinct 
from a state without that twist. In the ground state sec- 
tor without twists the states just have to be symmetrized 
in flavor; on the other hand, in sectors having twist loops 
the strings (quasiparticle tunneling loops) are altered, as 
demonstrated in the example of Fig. [9]:. 

Let us now give a physical derivation of the ground 
state degeneracy of the melted doubled Z 2 model on a 
torus. Before melting there were 4 2 ground states, la- 
beled by tunneling loops (non-contractible strings) of the 
quasiparticles: The 4 2 quasiparticles were just the un- 
constrained pairs (qp a ,QP/3) of Z 2 quasiparticles qp e 
{l,e,m, em} of the two flavors a, /3. In the melted 
phase there is a single ground state sector containing no 
twists. In this sector the quasiparticles must be sym- 
metrized in flavor, which gives us a total of 10 sym- 
metrized states: 4 states {qp a ,qPa) and 6 distinct states 
{qp ) + (qP0^ a ,qp a )- Next we consider a ground 

state sector containing one twist loop, say around the 
y direction. According to the argument depicted by 
Fig. [9j:, the twist effectively merges the two flavors, 
leaving us 4 ground states labeled by 4 Z 2 quasiparti- 
cles {l,e,m,em}. Finally, there are exactly 3 different 
twisted sectors, because the twist loop can wind around 
x, y, or both x and y (because there are only two fla- 
vors, higher winding numbers are irrelevant). All three 
twist sectors have the same degeneracy, because they are 
physically equivalent: The torus can be geometrically de- 
formed, without changing its topology, to map different 
twist loops onto each other (these transformations are 
Dehn twists, see Appendix [C]) . The torus GSD is there- 
fore 10 + 3 • 4 = 22. 

This result is easily generalized to the doubled Z„ the- 
ory on a torus. Namely, there are n 2 quasiparticles in a 
single copy of the model. The symmetrization in the un- 
twisted sector leads to n 2 + " (n 2 ~ ' ground states. There 
are still 3 twisted sectors on the torus, each contributing 
n 2 degenerate states. The total for the torus is therefore: 



\GSD me it2d,torus\ — TV 



, n A + 7 
'~2~' 



(25) 



confirming that for doubled Z 2 theory 
\GSD me it2d,torus(n = 2)| = 22. 



B. Physical argument for higher genus GSD 

Generalizing arguments of the previous subsection to 
an arbitrary surface of genus g is straightforward. Before 
melting, each of the g holes contributes non-contractible 
strings, leading to a total of (n 2 ) 9 degenerate ground 
states of a single copy of Z„ toric code. 

Each copy in the doubled model has n 2g degenerate 
states, which have to be symmetrized in the untwisted 
sector of the melted phase. This contributes 
degenerate states. 



Next, we consider the sector having a single twist loop, 
along the y direction of the g-th. hole (see Fig. [9]i) . Fo- 
cusing first on the states described by strings winding 
around and through the g-th hole only, we can apply 
the reasoning from the previous subsection. Namely, the 
twist merges the two copies of n 2 such states into only 
n 2 states which are automatically symmetric in flavor. 
The g-th hole therefore contributes a factor n 2 to the 
degeneracy of this sector. The strings tied to the other 
g — 1 holes contribute n 2 ^ 9-1 ) states per flavor, which 
have to be symmetrized. This gives a contribution of 

^ states from these holes. Combining all 

the holes, the total degeneracy in the sector with a single 

twist loop IS ' 7, J - J -. 

Again, because there are only two flavors, the twist 
loops that matter wind at most once. Each hole has two 
winding directions, giving the total of 2 2s possibilities for 
the winding of the twist loop. Exactly one of these pos- 
sibilities is when all the winding numbers are zero, i.e. 
the untwisted sector. There are therefore 2 2g — 1 twisted 
sectors, all of them physically equivalent, as argued pre- 
viously. 

Adding the untwisted and twisted sectors, the total 
GSD of the melted doubled Z„ toric code on a 2d mani- 
fold with genus g is: 



\GSD 



raeltld 



! = n 2g (n 2g + l) | (22fl ^(n^-D + l^ 



which reduces to the torus result in Eq. (25 1 when g = 1. 



C. The G n = (2„ xZ„) xi Z2 gauge theory: group 
structure, quasiparticle number and GSD 

We now turn to a gauge theory continuum descrip- 
tion of the melted lattice system. To start with, the 
doubled toric code is represented by the Abelian gauge 
group Z„ x Z„, a direct product of two flavor copies. As 
we argued in the previous subsections, the melting makes 
(|6| = a/2) dislocations dynamic, which means that also 
flavor-switching Volterra lines of all lengths occur in the 
system. This makes the two flavors physically locally 
indistinguishable. The flavor-mixing operation therefore 
becomes a local symmetry operation, i.e. a gauge trans- 
formation. 

According to these physical arguments, the gauge 
(structure) group of the melted doubled Z„ toric code 
system is formed by 



.92 








(27) 



where j a and jp label gauge transformations in the a 
and P flavored copy of Z„ gauge group, respectively; on 
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the other hand, the melting introduces the off-diagonal 
element gi which obviously exchanges the gauge fields of 
different flavors. The clement generates a Z 2 subgroup 
responsible for exchanging flavors, introducing a "twist" 
into the doubled Z„ gauge group we started from before 
melting. This twist in the gauge group is the mathemat- 
ical reason for appearance of "twist loops" we considered 
earlier (see Appendix [C|. The gauge group of the melted 
theory can be written as G n = (Z n x Z n ) x Z2, where 
the final Z2 acts by exchanging the first two factors and 
is generated by 52- Specializing to the n — 2 case, the al- 
gebra of above matrices explicitly gives the non-Abelian 
dihedral group D4 of order 8. 

The structure and physical properties of multi- 
component U(l) Chern — Simons gauge theory contain- 
ing a Z2 twist have been studied in Ref . 132] through the 
example of U(l) x U(l) x Z2. In the present case, be- 
cause the Z„ x Z„ gauge theory can also be described 
by a multi-component C/(l) Chern — Simons theory, one 
can apply the methods of Ref. S3] to study the proper- 
ties of the Z 2 twisted theory. In Appendix [C] we study 
properties of the melted phase, such as GSD, using this 
Chern — Simons approach. 

However, because the twisted theory in the present 
problem is described by a discrete gauge group G n , one 
can directly study the properties of the G n gauge theory 
More importantly, the multi-component U(l) Chern — 
Simons theory can only describe topological orders in 
2+ld. But the discrete gauge theory approach also ap- 
plies to higher dimensions, and we use this in Sec |VI| In 
the rest of this subsection we obtain the GSD in 2+ld di- 
rectly from considering the group structure of G n , and we 
confirm that result by solving the Chern — Simons theory 
in Appendix |C| 

The GSD of a topologically ordered state on arbitrary 
genus g surface, S g , can be determined directly from its 
quasiparticle content (for instance, see Refl43|): 



2(9-1) 



(28) 



where 7 labels the quasiparticles species, while d 7 is the 
quantum dimension of the 7 quasiparticle. On the other 
hand, when the topological order is described by a dis- 
crete gauge group G, such quasiparticle information fol- 
lows directly from the group structure of G. Namely, 
a quasiparticle is a dyon labeled by the pair 7 = (C, /i) , 
where C is a conjugacy class in G, physically representing 
the flux, and /x, representing the charge, labels an irre- 
ducible representation (IRREP) of the normalizer Nc of 
a representative element in CA- The quantum dimension 
of a dyon equals the product of class size (i.e. number 
of elements in C) and the dimension of the IRREP fi of 
N c - 

Table [TT] summarizes the relevant information about 
the group G n , which is easy to obtain using induction 
of representations from an Abelian subgroup. One can 
immediately establish from the table that there are in 



total: (a) 2n 2 quasiparticle species 7 having quantum di- 
mension dj = 1; (b) " *-" 2 — — species having dimension 
e? 7 = 2; and (c) 2n 2 species of dimension d 7 = n. Plug- 



ging into Eq. ( 28 ) , after trivial algebra it follows that the 



resulting expression precisely matches the GSD result, 
Eq. ( 26 ) , obtained from simpler physical reasoning: 



S g (G n ) = \GSD, n 



1 112<1\ 



(29) 



Since the genus g is arbitrary, it also follows that the 
quasiparticle content of the melted doubled Z„ theory 
matches the quasiparticle content of the G n topologically 
ordered theory. This is convincing proof of emergence 
of a non-Abelian G n topological order from a doubled 
Abelian Z„ order upon giving dynamics to Z 2 "flavor- 
mixing defects" (dislocations) as a general mechanism. 
In particular, the simplest melted doubled Z2 toric code 
is described by the non-Abelian D4 topological order. 



VI. EMERGENT NON-ABELIAN GAUGE 
THEORY IN DISLOCATION MELTED DOUBLED 
TORIC CODE IN 3D 

Even though the systematic study of three dimensional 
manifolds is an open problem, we can obtain some results 
by focusing on the simple case of a lattice with periodic 
boundary conditions, i.e. a model defined on the three- 
torus T 3 . 

Dislocation melting in 3d is a much subtler subject 
than in the 2d case, as we briefly touch upon in the Dis- 
cussion section. However, based on the simple physical 
insight that melting gauges the flavor symmetry in 2d, 
as described in the previous section, we conjecture that 
the melted phase of the doubled Z n toric code in 3d has 
topological order described by the G n = (Z„ x Z„) x Z 2 
gauge group. 

We test this conjecture only in the simplest case of 
n = 2, and as already mentioned, only on the three-torus. 
In this case, the melted phase of the Abelian doubled Z 2 
toric code is assumed to realize the non-Abelian G2 = 
D4 topological order. Using physical arguments as in 
Section |V A| we will derive the GSD of the melted phase 
to be 



\GSD melt3d (n = 2,T 3 



92. 



(30) 



On the other hand, in Appendix [P] we will prove that a 
D4 topologically ordered state on a three-torus T 3 has 
GSD also equal to: 



5 T s( J D 4 )=92, 



(31) 



thereby providing some evidence for the conjecture and 
the fact that in three spatial dimensions dislocation melt- 
ing of an Abelian phase can realize a non- Abelian phase. 

Let us now consider the doubled Z2 toric code in 
three dimensions with periodic boundary conditions, and 
find its GSD, which is \GSD meH3d {n = 2,T 3 )|. Before 
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No. of Cg 


C g labeled by g 


Class size, \C g \ 


Centralizer N g 


No. of IRREP in N g x IRREP dim 


No. of qp's 


n 


{x a y a }, o = 0...n-l 


1 


G n 


2n x ld + n(n 2 _1) x M 


n 2 (n + 3) 
2 


n(n — l) 
2 


{x a y b * a }, o,6 = 0...n-l 


2 


Z n x Z n 


n 2 x Id 


n J (n-l) 
2 


n 


{•wai"^}, a + i> = c; c = . . . n — 1 


n 


Z n x Z 2 


2n x Id 


2n 2 


Total number of quasiparticles in G n 


2 



TABLE II: Properties of (non-Abelian) gauge group G„ = (Z n x Z„) xi Z2 (n > 1), which describes the double-dislocation 
lattice melted phase of the doubled Z n toric code in 2d. An element g 6 G n is written as g = u s x a y b , where x, y generate the 
a, (3 flavored Z n , respectively, so that x n = y n — e and [x,y] — 0; the element u generates the flavor-mixing Z2, so ux = j/w, 
u 2 — e; therefore, s = 0, 1, and a, b = . . . n — 1. C fl is the conjugacy class labeled by a representative <? £ G n - The first 
column shows the total number of conjugacy classes of given type (3 types in total), and each class is labeled using parameters 
in second column. Centralizer N g is the subgroup of elements that commute with g. Every conjugacy class labels a magnetic 
flux. Quasiparticles are dyons, carrying a flux and a charge, where the charge is an irreducible representation (IRREP) of N g , 
with Cg labeling the flux. Hence, column six shows the total number of dyons carrying flux of the type labeled by the table 
row; it equals the product of the first column and the first factors in fifth column. The quantum dimension of a dyon equals 
the product of C g class size (third column) and the dimension of the IRREP of N g (second factors in fifth column). 



melting, each flavor has GSD equal to 2 3 , due to non- 
contractible strings. 

The melting enforces symmetrization of flavors, and 
also introduces "twists". In 3d, the twist occurs across 
a plane that spans the system. The twist plane can be 
imagined as created by a dislocation loop which is nucle- 
ated and stretched until it annihilates using the periodic- 
ity of the system. A non-contractible string on the three- 
torus has to switch flavors as it pierces the twist plane and 
spans the periodic system (see the horizontal string in 
Fig. [9J: for a 2d analogy). Since there are only two flavors, 
only strings which span three-torus directions X, Y, Z up 
to once are relevant. There are therefore 7 twisted sec- 
tors, each having a flavor-switching string which spans 
X, Y, Z, XY, YZ, XZ, or XYZ. The twisted sectors 
are physically equivalent. Analogously to the 2d case, a 
twisted sector has the two flavors replaced by a single 
Z2 copy, which is automatically flavor-symmetric, and 
therefore has a 2 3 degeneracy. The twisted sectors there- 
fore contribute 7 ■ 8 = 56 states. The untwisted sector 
has to be flavor-symmetrized, therefore 2 3 x 2 3 leads to 

36 states. The total GSD is 92, as claimed 



a) 



2 J (2 J ~1) 



in Eq. (30 1. 



VII. DISCUSSION AND CONCLUSIONS 

In this paper we explore the dependence of ground 
state degeneracy in simple exactly solvable topological 
models (in both 2d and 3d) on the complicated topology 
of the underlying space manifold, where that topology 
is provided by dislocation points (in 2d) and arbitrarily 
linked and knotted dislocation lines (in 3d) with Burg- 
ers vector connecting two translated copies of the model. 
Because we are studying topological properties, our re- 
sults are valid for the entire phases to which the exactly 
solvable models belong. 

The main new physics in this paper is about the 3d 
models. In 3d, if the ground state degenearcy(GSD) de- 
pends on the knotting/linking of dislocation loops, the 



d) 



c) 



FIG. 10: Three-dimensional toric code in a system with open 
boundary conditions. Dislocation lines (blue) and flux strings 
(black solid/wavy for two copies) can end at the surface, 
(a) The string loop between two dislocations would be non- 
contractible if it were not for (b) the surface and the fact that 
(c) string endpoints on the surface can avoid dislocation line 
endpoints. Any number of open dislocation lines beyond the 
first cannot contribute to topological GSD. (e) String loop be- 
tween dislocation loop (in bulk) and open dislocation line is 
a non-contractible string. The open string can never be con- 
tracted because it (f) intersects the dislocation loop Volterra 
surface an odd number of times. 



GSD can only be a topological invariant of the links or 
knots of the dislocation loops. However, we find that the 
doubled Kitaev's toric code shows no dependence of the 
GSD on the knotting properties of dislocation loops. In 
this 3d topologically ordered phase, GSD only depends 
on the total number of (possibly linked) dislocation loops. 

In fact, it turns out that in the limit of large number 
k of dislocation loops (pairs in 2d), the GSD scales as 
2 k in 3d. As a dislocation loop can always be allowed to 
be added into the lattice, the extra GSD introduced by 
one dislocation loop should be an integer. In this context 
therefore the doubled Kitaev toric code is reasonably the 
minimal theory in 3d in which dislocations introduce ex- 
tra GSD. In 2d however the scaling is 4 fc , indicating that 
there might be a more fundamental 2d theory, a "square 
root" of the doubled 2d Kitaev toric code, in which the 
GSD scales as 2 k . Such a "square root" theory is indeed 
realized in the plaquette model version of the 2d toric 
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cod<P, and 2 fe scaling of the GSD was found. 

T he Ki taev toric code can be generalized to Z„ toric 
cod d 45 * 46 -!, whose low energy effective theory is described 
by the Z„ gauge theory, and when n = 2 the Kitaev toric 
code is restored. Our results can be generalized without 
too much difficulty to doubled Z„ toric code in both 2d 
and 3d. Here we present the main results and sketch the 
strategy of proof. The detailed analysis and proof can be 
found in Appendix |A"| 

Obviously, in 2d in a doubled Z„ toric code with pe- 
riodic boundary conditions with k pairs of dislocations, 
\GSD\ = n 4 ■ n 2k ~ 2 which is a direct result of the ge- 
ometrical interpretation of dislocations in 2D (see Fig{2] 
or Ref. 12). A "square root" of such a theory has been 
recently studied in RefHTl where GSD scales as n k when 
k is large. 

In 3d, we find that the GSD in a doubled Z n toric 
code with periodic boundary conditions in the presence 
of k dislocation loops is \GSD\ = n 6 ■ n k ~ 1 , independent 
of linking and/or knotting of the dislocation loops. The 
proof of this result is a slightly non-trivial generalization 
of our proof for the doubled 3d Kitaev toric code. Fol- 
lowing our strategy of using string/membrane operators 
to resolve the ground state degeneracy, one still can con- 
struct the string operators which can be interpreted as 
"electric flux lines " in a 3d Z 

n S&rige theory. But in a Z n 
gauge theory with n > 2, the flux line has a Z„ valued 
strength. If we focus on the flux line with the fundamen- 
tal unit of strength (strength=l), then the flux line has 
a well-defined orientation. Some essential steps in our 
proof for the doubled Kitaev toric code now need to be 
reconsidered carefully. For instance the step 6 illustrated 
in Fig |6] requires annhilation of two string operators. But 
in a Z„ theory two string operators both with strength=l 
cannot annihilate if n > 2. In Appendix [X] we carefully 
reconsider the rules of surgeries involving the string oper- 
ators and the proof still goes through (basically because 
dislocations allow a string of strength s to be converted 
to strength n — s, which is its inverse). 

Can the ground state degeneracy depend on linking or 
knotting of the dislocation loops in a 3D theory then? 
Our study of the Z n Abelian gauge models strongly sug- 
gests that at least for Abelian gauge theories this may 
be impossible. However for non-Abelian gauge models, 
for instance the discrete non-Abelian gauge theories with 
deconfined phases in 3D, our proof based on string op- 
erators may be invalidated and a knotting or linking de- 
pendent GSD may occur. We leave this interesting pos- 
sibility as a subject of future investigation. 

What happens in a system with open boundaries? It is 
well known that the 2d Z2 toric code can have two types 
of edges: An open electric (magnetic) flux line ending on 
first (second) type of edge commutes with the Hamilto- 
nian, while the open magnetic (electric) can not commute. 
The types of edges in the system therefore control the 
number of non-contractible electric and magnetic lines 
(strings). (These edge types are realized by different ter- 
mination of the lattice.) However, it does not make sense 



to put dislocations at the edge of the 2d system, and it 
is easy to see that the dislocations add the same extra 
GSD in open as in periodic 2d systems. 

In 3d, one can similarly terminate the lattice model 
at the surface of the system in two ways: One allows 
open strings to end at the surface while still commuting 
with the Hamiltonian, while membranes terminating at 
the surface cost energy. The opposite is true for the other 
surface type. Let us focus on the first type, for simplicity. 
Firstly, any non-contractible string piercing a Volterra 
surface of a dislocation loop in the bulk will stay non- 
contractible even when the system has a surface, since the 
number of piercings of that Volterra surface will remain 
odd (and the string therefore non-contractible) no matter 
how we manipulate the string using the surface of the 
system (see example in Fig. 10 2,f ) . 



The situation is different for an open dislocation line 
that ends at the surface of the system (it has two end- 
points on the surface). As Fig. 10 1 demonstrates, the 



string loop that pierces Volterra surfaces of two such 
open dislocation lines becomes contractible. Namely, the 
string loop can break into two open strings by reaching 
the surface, and the string endpoints can avoid the dislo- 
cation lines' endpoints on the surface, allowing the open 
strings to contract (Fig.|l0^-d). This immediately shows 
that any (non-zero) number of open dislocation lines is 
equivalent to a single one, as far as extra topological GSD 



is concerned. However, as Fig. 10 5,f clearly demonstrates, 
the argument in the previous paragraph ensures that a 
single open dislocation line still behaves as a closed dis- 
location loop in the bulk, because a string that pierces 
through Volterra surfaces of a dislocation loop and an 
open dislocation is still non-contractible. 

Therefore, in a 3d system with a boundary allowing 
commuting open strings, the extra GSD due to k dislo- 
cation loops in the bulk (i.e. 2 fe ~ 1 ) is multiplied by a 
factor 2 if there is any (however many) open dislocation 
lines that terminate on the surface. 

The extra topological GSD due to dislocations leads 
to additional zero-temperature entropy. This entropy is 
obviously extensive in the number of dislocation loops 
in the bulk, and therefore could be of experimental rel- 
evance. Let us emphasize that the GSD properties hold 
for the entire topologically ordered quantum phase, be- 
yond the exactly solvable doubled model. As sketched 
above, the extra entropy should stay extensive in num- 
ber of bulk dislocations even for realistic samples with 
surface. 

Recently, topologically ordered phases which also have 
global symmetries, i.e. symmetry enriched topological 
(SET) phases have been studied intensively. The topo- 
logical order is usually described by gauge symmetry and 
understanding its interplay with the global symmetry is 
the central issue in studying SET phases. In this work, 
non-Abelian topological order arises due to the promo- 
tion of the global Z2 symmetry (the flavor symmetry) to 
a gauge symmetry. Thus defect condensation seemingly 
provides a physical way to make the connection between 
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global and local symmetry, and therefore between differ- 
ent SET phases. 
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Appendix A: Generalization to Z n toric code 

Here we present the generalization of our GSD proof 
based on string operators of the 3dTC (a Z 2 theory) to 
a Z n version of the toric code. The Z n generalized toric 
code was introduced in Refl4"51 and we use a description 
of it similar to the 2d case from Ref. |46] In this way 
the lattice and stabilizer definitions are similar to the 
3d Z2 toric code of the main text. We first introduce 
the building blocks and then identify the strategy for 
extending our GSD proof using string operators. 

At each spin site i of the cubic lattice (see Fig. [3]) of 
the Z„ model we introduce n-state degrees of freedom 
\qt) , which are acted upon by generalizations of the Pauli 
matrices of the Z2 model. These unitary operators Xi, Zi 
at site i form a quantum rotor algebra: 

Zi = LO q ' I Qi) , 

Xi \<li) = \<li - 1 mod n) , 



x\x t = z\Zi = 1, 
XiZi = LoZiXi =>■ XiZ\ = uj* ZjXi 
to = e l2v/n , 



(Al) 



with to* the complex conjugate. These relations general- 
ize the anticommutation of Pauli matrices. 

The stabilizers of the model of flavor a are presented 
in Fig. 11 and are defined similar to Eqs. ( 15|16 ): 



^4"g{ } = A"' " — X s + x s + ij /2X s + z /2X s-X /2X s -y ^X g_ 

(A2) 
(A3) 



au Ae,a _ ao, at 

s£{e} — A s — A s > 



B" 



Z 



z. 



>z 



(ab)+a/2 p(ab)+b/2 p (ab)-a/2 p(ab)-b/2> 



(A4) 



with the key difference being the necessary bipartite de- 
composition of the cubic lattice such that even ( "e" ) and 
odd ("o") cubic sites have star operators Af a and A°' a , 
respectively, positioned on them. This property also oc- 
curs in the 2d Z ra modePSl. Since the even and odd star 
stabilizers act with inverse operators on the spins they 




FIG. 11: Three-dimensional doubled Z n generalization of 
toric code. See caption of [3] for the basic ingredients shared 
with Z2 toric code. Three /3 flavored plaquette operators are 
shown, acting with Z (blue wavy lines) and Z^ (green wavy 
lines) operators. The operators are chosen so that the lo- 
cal cubic constraint for the plaquettes holds. Star operators 
on the "odd" cubic sublattice act with X operators (blue zag 
lines), and with X^ if on the "even" sublattice (only four stars 
of /3 flavor, positioned in one plane of the lattice, are shown) . 



share (Xi are unitary), the global constraint 

n m=* 

sG{e,o} 



(A5) 



remains as in Eq. (J5J). A plaquette stabilizer consists of 
operators Zi which alternate between daggered and nor- 
mal as one goes around the face of a cube which defines 
the plaquette, so that opposite sides of the plaquette both 
carry either 2 or 2*. This ensures that any plaquette 
shares two spin sites with a neighboring star such that 
both Zi and zj appear, while the star contains all X 
or all X* operators, leading to the commutation of the 
s/2plaquette and star according to Eq. (Al). The choice of 
which sides of a plaquette carry the daggered operators 
is made in Eq. ( A4 ) such that local cubic constraints are 



satisfied as before (see text after Eq. (|16|)). 

The j3 flavor of the model is created by superimposing 
a copy of the lattice translated by y/2 (Fig. [3]), as before, 
and the Hamiltonian becomes 



f,p 



v 



v 



(A6) 



Because of the commutation and unitarity, it is enough 
to count the stabilizers A, B and fix their eigenvalues in 
the Z„ theory. Since the constraints carry through as 
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string loop 
start point (not 
commuting) 



FIG. 12: String operators in Z n doubled toric code with dislocations. Only a single plane of 3d lattice is shown for clarity. The 
stabilizers are naturally repaired, because the dislocation cannot directly mix the even/odd sublattices; example shown: even 
and odd star. Even/odd sublattice is mixed after encircling the dislocation twice: the string constructed step-by-step clockwise 
on the right of the figure, starting from the star labeled by magenta is geometrically closed (flavor labeled by full/dashed line is 
switched twice), but does not commute with all stabilizers! Namely, the string has alternating Z ' j 'Z' operators (slanted black 
lines of different angle) to ensure commuting with stars, but this fails for the star at the start point. On the left of the figure 
is an example of a valid (c omm uting) string loop which changes flavor and also gets inverted, as demonstrated by its arrows 
(see paragraphs after Eq. (A8l). This loop is contractible. There is no inversion when the flavors on the top/bottom half of 
the loop are chosen oppositely, as depicted in Fig.[l3p(left). 



for the Z2 toric code in the main text, the GSD counting 
remains the same, except that each unconstrained "spin" 
degree of freedom now carries n-fold degeneracy. This 
gives \GSD 



(n) 
3dAdeal I 



Our next step is to introduce string operators consis- 
tently in the Z„ theory with dislocations. A string oper- 
ator is non-trivially different from the Z2 case, since its 
eigenvalue, i.e. the flux it carries, can be any of n values 
e" 27T /", s € {0, . . . , re—l}. The string must be defined us- 
ing both daggered and normal operators: given a closed 
path (loop) P/ which stepwise connects neighboring star 
sites of flavor /, the Z operator acting on the spin be- 
tween two stars has to alternate between daggered and 
normal with each step: 



-i-r \Zi, i odd, 

itf f W' ieven ' 



(A7) 



compare to Eq. (22). This way, a star operator crossed 
by the path Ft will share exactly a Z and a Z' operator 
with the string, therefore commuting with it. 

The crucial quality of the Z„ theory is that al- 
though the string Y,f(¥f) carries some flux e Js27r /™, s G 
{0, . . . , n — 1}, i.e. it has strength s, it is easy to define 
its inverse: 



E'(P / )=E'CP/) t , 



(A8) 



which of course carries flux n — s while having the exact 
same geometrical shape. The string is canceled by its 
inverse, but not necessarily by itself, in contrast to the 
case of Z2 theory (see Fig. |l3|a)). 



Our main proof for the GSD in presence of dislocations 
crucially depends on canceling two copies of a string of 
a given flavor (see step 6 of Fig. |6|. In the Z„ theory, 
starting from a string of strength s, the two copies of the 
string would not necessarily cancel. 

However, the proof is rescued by a remarkable fact: the 
aj3 surface of the dislocation allows inversion of a string, 
as well as its flavor change (Fig. \l3\p ). As we will next 
show, this remarkable fact hinges on the dislocation dis- 
rupting the even/odd division of the cubic lattice, which 
itself is directly related to mapping a string to its inverse. 

It is important to notice that the considered disloca- 
tions with |6| = a/2 directly couple the flavors, but do 
not simply couple the even/odd sublattices which are ac- 
tually connected by an a translation. The fact that the 
dislocation can not directly exchange even and odd sites 
also ensures that the stabilizers can be smoothly and con- 
sistently "surgically" repaired in a natural way across the 
a/3 surface (see next paragraph and Fig. [l2]) . 

Let us now show the crucial properties of a string 
and its inverse in relation to dislocations (Figs. 12|13 1. 
For readability, figure [12] shows only a single layer of 
the 3d lattice (which is similar to a 2d model), but the 
conclusions hold generally in 3d. As is obvious from 
Eqs. (A7l,(|A8j), the inverse of a string is obtained by 



switching the odd and even steps in the string opera- 
tor. The removal of star sites by the dislocation's af3 
surface can play exactly this role, as is explicitly shown 
in Fig. |12| Given a spin site i on the step between an 
odd and even star site (which are always well defined, as 
the spins they act on must be X xor X 1 * even after lo- 
cal surgery due to dislocation) , the string and anti-string 
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a) 



b) 






O 



FIG. 13: Continuum rules for strings in the 3d Z n dou- 
bled toric code, (a) A string carries a flux value labeled 
by s £ {0, — 1}. In presence of dislocations, an ar- 

row is a useful local representation of the flux carried by the 
string, because the string can get inverted (s — > n — s). In- 
verse strings cancel, while doubling a string only doubles its 
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strength (s — > 2s). (b) Top view shown for clarity, see Fig. 
for explicit lattice demonstration. (Left) Dislocation a/3 sur 
face changes flavor, but also can invert the string (s — > n — s). 
(Right) Only winding w — 4k strings commute with stabi- 
lizers (here k = 1 is shown), since a |46| = 2a translation is 
needed to avoid a jump between even/odd sublattices. These 
winding loops are the generalization of w — 2k ones from the 
Z2 model, and are also naturally locally contractible. (c) Us- 
ing rules from part (b), it is easy to show that dislocations can 
locally independently invert or change flavor of string loops. 
This suffices to ensure the possibility of step 6 of the GSD 
proof (see Fig. [6| within the Z n theory, upholding the entire 
construction. 



differ by having Z\ or Zj. Therefore, we are able to con- 



up to a translation by \2b\ — a, which is a jump from 
even to odd and vice versa in the division of the cubic 
lattice (for both flavors). One expects that it is required 
therefore to encircle the single dislocation line four times 
to obtain a valid string loop. This is indeed true. Fig- 
ure [12] explicitly shows how encircling a dislocation line 
twice, and therefore piercing the a/3 surface twice in the 
same direction, produces a geometrically closed loop (by 
switching flavor twice), but also gives a non-commuting 
string! Namely, the star stabilizer (being all X or all 
X^) at the start /end point of the path shares two sites 
with the string, both on which the string is Z (or Z'). 
(See commutation relations in Eq. ( Al ).) The successful 
generalization of the double-winding string loop to the 
Z n theory gives the 4-winding version; this kind of loop 
is also naturally expected to be contractible using the 
locally surgically repaired stabilizers at the dislocation 
line, and is shown in Fig. |13p . 

Using only the properties of string loops explained 
above, one can find the remarkable (and useful) gener- 
alization of rules from Fig. [5] The crucial rule is that a 
string encircling locally parallel dislocation line segments 
with an a/3 surface stretched between them can, by using 
only local operations, change its arrow direction, as well 
as its flavor. This rule, in the continuum representation, 
is depicted in Fig. 13 :, and is the generalization of rule 
in Fig. 5b. Concerning other rules, a proof of the one in 



Fig. 13 D(left) is shown explicitly on the lattice in Fig. 12 
These rules allow one to use all the constructive steps 
in our proof (Fig. [6]), with slight modifications. Firstly, 
the flux s has to be assigned to the string loop. Once 
the phase between steps 3 and 4 is reached, one must 
also explicitly state the arrow direction for the loop, that 
accompanies its strength s. It is easy to see that the two 
copies of the string in step 5 are not only geometrical 
copies, but they also have the same arrow direction, as 
well as the same strength s. It is simple then to use 
the new string rules to achieve the inversion of one of the 
copies (i.e. switching it to strength n—s). The two copies 



struct a simple pictorial way to locally distinguish the then cancel, just as is required by the original proof. 



string and inverse string, exemplified in Fig. 12 We can 
assign an arrow at each spin site along the geometrical 
line following the string, i.e.: following the path of the 
string loop in a fixed orientation (call it "clockwise"), if 
on the step between odd and even star site the string acts 
with Z operator, assign an arrow pointing along the path 
(so in the "clockwise" direction). Otherwise assign the 
opposite direction. One can easily see (e.g. using Fig. 12 ) 
that the even/odd definition of stars and strings will pre- 
serve the arrow direction all along the string, at least in 
absence of dislocations. Also, deformations of the string 
by the action of plaquettes (and/or anti-plaquettes B^) 
must preserve the arrows due to plaquettes' alternating 
Z/Z^ geometry. By this construction, two loop segments 
sharing a geometrical path but having opposite arrow 
directions on them cancel each other, and therefore rep- 
resent the string and anti-string on that path segment. 

Intuitively, encircling the dislocation line twice adds 



Appendix B: The membrane operator in presence of 
dislocations 



Here we complete our analysis based on string oper- 
ators. Strictly, a non-contractible string has to be ac- 
companied by a non-contractible membrane it pierces an 
odd number of times, for our claims about the GSD in 
subsection IIV Dl to be true. 

We first consider the practical construction of mem- 
branes on the lattice, and introduce the useful concept of 
"minimal membrane" . Then we present the analysis in 
the continuum. 

A particularly important membrane is one that con- 
tains the entire chosen dislocation loop. As we check 
directly on the lattice (see Appendix [E| using star oper- 
ators, and discuss further later, such a membrane cannot 
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FIG. 14: Elementary operations on membrane operators in presence of dislocation loops, derived from their behavior on lattice. 
Blue lines are dislocation lines, and orange straight/wavy surface is a membrane operator of a//3 flavor, (a) The minimal, 
locally non-contractible membrane (red line shading) has two halves of opposite flavor, tightly spanning a dislocation loop. It 
can expand into membrane of either flavor, (b) Linking dislocation loops preserves minimal membranes; local redefinition of 
two minimal membranes is only needed where two a/3 surfaces intersect (green line), see text, (c) Local passage of a piece of 
membrane through dislocation loop, (d) Relevant minimal membranes do not have to coincide with a/3 surfaces. An example 
of two dislocation loops is shown, with different choices for the a/3 surfaces. The non-contractible strings that anticommute 
with the chosen minimal membranes are also shown. 



be contracted locally, i.e. through the dislocation loop 
it contains. (In absence of other dislocations, the mem- 
brane can be expanded until it is contracted to nothing 
due to the periodicity of the lattice — this is a "global" 
operation, as performing it depends on whether there are 
other dislocations in the lattice.) Using the properties of 
the enclosed a/3 surface, we show that the membrane in 
fact can be locally transformed to the opposite flavor. 
It is numerically favorable to have a simple prescription 
for such a membrane (whether actually it is contractible 
or not in given circumstances) such that it is locally of 
the smallest possible volume. Remarkably, such "mini- 
mal membranes" are easily constructed even in the case 
of linked dislocation loops (when two a/3 surfaces inter- 
sect, Fig. 20 1. The prescription for constructing minimal 



membranes in presence of various dislocation loops fol- 
lows: 

(Edge loop) Apply all star operators that are going to 
be removed in creation of the a/3 surface; after site 
removal, the leftover a x operators form the mini- 
mal membrane; these are actually acting on nearest 
neighbor (opposite flavor) spins on the two sides of 
the a/3 surface. Summarily, 



- edqemin 



TIa. 



(Bl) 



where V is the set of star operators that are re- 
moved by the Volterra (a/3) surface. 

(Screw loop) Apply all star operators that act on spins 
in screw plane, then apply the new, repaired ver- 



sion of same stars after screw operation; all these 
a x operators form the minimal membrane; these 
are actually acting on all spins in the screw plane 
(both flavors), except along the dislocation edge- 
B (Fig. 19 1 where the new stars are more heavily 



changed by local surgery (Fig. 19). The procedure 
can be summarized by the formula: 



scre wm t ii 



n a * a *> 



(B2) 



s£V 



where the set V contains all star operators that act 
on spins in the Volterra surface, and A s are stars 
in the original lattice, while A s are the same stars 
after "surgery". 

(Linked edge and screw loops) Apply screw loop 
surgery of lattice; construct the screw minimal 
membrane (will be edited in last step); construct 
edge loop minimal membrane (note that in this step 
the stars on the intersection of the two a/3 surfaces 
have already been edited by screw loop); finally, 
apply edge loop surgery on lattice (which removes 
some operators from the screw membrane) , Fig. 20 



The last case constructs two minimal membranes, for 
each of the linked dislocation loops. These two mem- 
branes are actually independent, i.e. they cannot cancel 
each other. In other words, however we modify the prod- 
uct of these membranes by multiplying it with stars in 
the final dislocated lattice, there will always remain a x 
operators acting on some spins, and thus the two mem- 
branes cannot cancel each other. 
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To analyze the membranes as we did the strings, we 
have directly inspected all situations on the 3dTC lat- 
tice (Figs. 16p0 ) to derive a set of general rules that 



the membranes obey as they are deformed and moved 
through the lattice by the action of stars which have un- 
dergone proper "surgery" due to the dislocated lattice 
(Appen dix [E| . Having in mind the detailed discussion in 
Section |III A| we fix the Volterra a/3 surface of each dislo- 
cation loop, and retain the global flavor assignment from 
the original ideal lattice. 

The basic set of rules for the behavior of membranes 



is presented in Fig. 14 which also shows their relation to 
non-contractible strings. Intuitively, a membrane is an 
elastic closed surface which can pass through dislocation 
loops. However, if it is created such that it encompasses 
an entire dislocation loop (and its a/3 surface), it cannot 
locally be contracted, but instead can just locally (using 
the a/3 surface) change flavor. Such a membrane we will 
call "minimal" (for the given dislocation loop). 

Following our analysis for strings above, we again con- 
sider the two simple examples: 

• A single dislocation loop does not change the GSD: 
The minimal membrane is trivial, as it can be ex- 
panded and, using system periodicity, contracted 
far away from the dislocation loop, as no other dis- 
location loop is in the way. 

• Two separate dislocation loops add one non-trivial 
membrane. Here the minimal membrane of one 
loop can be expanded in the periodic system un- 
til it engulfs the other loop tightly. So, the two 
separate minimal membranes of the two disloca- 
tion loops are actually equivalent. Further, due to 
the flavor inversion rule the two flavors of the mini- 
mal membrane are also equivalent. We are thus left 
with only one independent non-contractible mem- 
brane candidate. The standard system-spanning 
strings must pierce this closed membrane an even 
number of times (possibly times). The new non- 
contractible string (which pierces both a/3 surfaces 
once) is topologically forced to anticommute with 
the minimal membrane, because the membrane can 
be considered as of single flavor (either one) , while 
the string has a flavor change on the a/3 surface 
which is inside the membrane. Notice that using 
system periodicity to change the meaning of word 
"inside" to "outside" gets us nowhere, as this oper- 
ation topologically only maps one dislocation loop 
onto the other. 



minimal membranes. However, this does not mean it is 
dependent! Consider the old system, before adding the 
new dislocation loop. A product of all minimal mem- 
branes there is globally contractible, by expansion away 
from the dislocation loops, and using the periodicity of 
the 3d system. With the new loop however, expanding 
the product of all old minimal membranes does not con- 
tract it, because the new loop is in the way. We just 
establish that the new minimal membrane is equivalent 
to the product of the old ones, but this is not a triv- 
ial membrane anymore. Therefore the product of all the 
old minimal membranes is now a non-contractible oper- 
ator that anticommutes with the new string, due to the 
string's single piercing of the new dislocation loop. 

From the membrane perspective it is even more re- 
markable how dislocation line linking does not alter the 
GSD. The minimal membranes of two dislocation loops 
stay well defined and distinguishable operators, even as 
the two dislocation loops are pushed into each other, and 
their a/3 surfaces intersect (see Fig. 14 3, and the ex- 



plicit lattice construction further in this section, as well 
as Fig. 



20). The two linked minimal membranes stay 



separate in the following sense: If we start from minimal 
membrane of the first a/3 surface, expand it using al- 
lowed deformations by stars, then once it grows through 
the periodic system, it will shrink back precisely into the 
form of the minimal membrane for the second a/3 sur- 
face, which belongs to the other linked dislocation loop 
(assuming no other dislocation loops in the system). 

Finally, let us clarify the 3dTC with parallel disloca- 
tion pairs spanning the system in one periodic dimen- 
sion (say Z). In this case the membranes orthogonal 
to the XY plane take the role of the 2dTC a x strings, 
due to translational invariance along z. The counting in 
the 2dTC therefore matches the 3dTC when we consider 
any XY cut through the 3d system. Additionally, the 
3dTC still has a single Z spanning string (two flavors are 
equivalent) and a single XY membrane (the two flavors 
are equivalent), so that there is one extra anticommuting 
pair in the 3dTC, leading to the 



\GSD. 



3d, spanning \ 



2\GSD 2d \ 



(B3) 



See Table U for exact numerical results that corroborate 
this, and Section |IVD| for additional details about the 
strings. 



Moving on to k dislocation loops, one needs to consider 
an addition of one new dislocation loop. The appropri- 
ate candidate for the new membrane must be the mini- 
mal membrane of the new loop. Notice how the situation 
here is different than for strings. Namely, it is immedi- 
ately obvious that by expanding in the periodic system, 
the new membrane just encompasses all old loops, and 
is therefore just the combination (product) of all the old 



Appendix C: Computing the ground state 
degeneracy of G„ gauge theory using 
multi-component (7(1) Chern — Simons theory with a 

twist 

The G n = (Z n x Z„ ) xi Z 2 gauge theory can be de- 
scribed by the multi-component U(l) Chern — Simons 
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theory with the following K-matrix: 



K 



K a =K fi = 




(CI) 
(C2) 



with two blocks for two flavors (copies). The Chern — 
Simons theory on the space manifold M has the La- 



grangian 



n 
An 

1 

47T 



M 



f K IJ A I dA J 
Jm 



(C3) 



(C4) 



a 4-vector con- 



with A = A n dx», A„ = (a«, a°, < aft) 1 
sisting of two flavored parts with components 

A{ = {al~ai) T , (C5) 

while x^ — t,x,y. The Lagrangian needs to be supplied 
with proper boundary conditions; crucially, the gauge 
theory contains flavor-exchanging Z 2 twists. On a torus 
the twists are represented by 

A fl (f+ L&) = (of <8 \i)A il {r){u\ i <g> 1 2 ), (C6) 

where e Xl e y = 0, 1 determine whether there is a Z 2 twist 
along the given direction, and L, is the system size in 
direction i = x,y. The theory on the torus therefore sep- 
arates into 4 sectors (one untwisted, three with a twist). 
On a genus g surface, the presence of twists along the 
2g elementary cycles (non-contractible loops) selects one 
out of the 2 2g topological twist sectors. 

In the next two subsections we evaluate the ground 
state degeneracy on the torus (\GSD\ = and a 



genus g surface (\GSD\ = 
the nature of the theories. 



(n) 



as well as comment on 



Let us note in advance that the Sg result for G n de- 



rived below, Eq. ( |C37| ) (or Eq. (|29| and Eq. (|26j) of main 
text), formally coincides with the result for twisted (sin- 
gle flavor) Z„2 theory described by Z„2 x Z 2 Before the 

twist the latter theory has a 2x2 if-matrix, K = ( ® 2 'q 

(contrast to Eq. (Cf )). This implies that the number of 
species and quantum dimensions of quasiparticles also 
coincide between G n and Z„2 x Z 2 , remarkably since the 
theories are fundamentally different, and the difference is 
manifest in the derivation of Sg itself (compare Ref. 02] 



to section C2|. However, due to the physical argument 
of the gauge structure (see Sec V C I , we expect that the 
G n is the correct description for the physical state of the 
melted doubled toric code. 



1. Degeneracy on a torus 

We first consider the ground state degeneracy of the 
(Z„ x Z„) x Z 2 theory of melted toric code on a torus 



{L, 



L for simplicity). Setting the gauge Aq = 



implies that the gauge fields of flat connections are irro- 
tational and characterized by their holonomies along the 
non-contractible loops on the torus. With an additional 
rescaling, 

Af(trf = ^{Xt(t),Xf(t)) T , (C7) 

Af(trt = j-{Yf(t),Yt{t)) T , (C8) 

these uniform fields therefore parametrize the gauge- 
inequivalent configurations. We additionally need to take 
into account that large gauge transformations enforce the 
identifications 



(X f + k f ,X f + k f ), (C9) 
(Y f + lf,Y f + If), k f ,i f ez. (CIO) 



(x f ,xf) 

(Yf,Y f ) 

The Lagrangian from Eq. (|C3j) becomes 

L = 27m^ (x f Y f +X f Y f ^ , 
/ 



(Cll) 



leading to a vanishing Hamiltonian and therefore a wave- 
function that is a linear combination of simple plane 
waves: 



il){X f ,X f ) = ^2 C^ZY[S[2TrnX f -2np f ]S[2T:nX f ~2Trq f l 

p a ,q a f 

(C12) 

with pf, qf being integers, pf, qf = . . . n — 1, and S[x] 
the Kronecker delta function. We can label states by 



\ip) = \p a ,q a ; 



(C13) 



and gauge equivalence from Eq. ([C9| limits the number 
of independent coefficients which enumerate the indepen- 
dent states. Namely, 



p a —n,q a 

^p a ,q a 



(~iPp,W 

Pa ,qa~n 

(-iPftm- 11 



(C14) 



giving n 2 independent states per flavor, since Pf,qf — 
0...U-1. 

The independent states in the untwisted sector are 
physical only if they are Z 2 symmetric. This constraint 

leads to n 2 states of the form \p,q;p,q) and " ^ 
states of the form \p a , q a ; pp^qp) + \pp, qp;p a ,q a ), for the 
total of 



g=l , untwisted 



n 2 (n 2 + I) 
2~ 



(C15) 



There are three twisted sectors, having a Z 2 twist along 
x, y, or both directions of the torus. There exist symme- 
try transformations mapping the sectors onto each other, 
so that all three have the same degeneracjEl. Focusing 



2G 



on a y-twist example, in accordance with Eq. ( C6 1 we 
demand that 



Af(x,y + L)=A^(x,y) 
A?(x,y + L)=A?(x,y), 
A{{x + L,y)=A{(x,y). 



(C16) 
(C17) 
(C18) 



We can now switch to a field B? defined on the system 
with double length along y. Physically, we glue the two 
tori of opposite flavors together after cutting them open 
completely along one direction (y) , which is along the dis- 
location "tunneling loop" winding along y. (This loop is 
the remnant a/3 line of the annihilated dislocation — anti- 
dislocation pair.) The resulting system is topologically 
still a torus. The new field is constructed as: 



B n(x,y) = {b fl ,b fi ) T = 



Afe,y-L), 



< y < L 

L <y <2L 
(C19) 



and satisfies 

Bpfay) = B M (x + L,y) = B M (x,y + 2L). 

Similarly to the untwisted case, we can set Bq 
parametrize 



Bi(t,x,y) 



2tt 



(X i (t),X i (t)) T , 



(C20) 
and 

(C21) 



where L y — 2L X = 2L. The Lagrangian Eq. (C3| be- 
comes 



n 
An 



dx \ 

ii Jo 



2nn X X X 2 



dy (bdb 

X1X2 



bdb) = 



(C22) 



while the large gauge transformations consistent with 
Eq. (C19) introduce the periodicity 



Xi 



Xi + 1. 



(C23) 



The trivial equations of motion and periodicity lead to 
the wavefunction 

jj,(Xi) =^2c Ptg 5[2imX 1 -2wp]5[-2nnX i -2'irq], (C24) 

where c Piq = c p ~ n . q = c Piq + n , and therefore the degener- 
acy of the y-twisted sector is n 2 . Note that the twisted 
sector is equivalent to a single flavored (Z 2 invariant) un- 
twisted sector (e.g. Eqs. ( |C22| ) and ( |C11[ >). 

Since there are 3 twisted sectors, it follows 



S 



(n) 

g—l,twisted 



3n 2 , and the total ground state degen- 



eracy on the torus is 



rr(n) 

- 



n 2 (n 2 + 7) 



This result precisely matches the one obtained from a 



simpler physical argument for the melted state, Eq. (25 1 



2. Degeneracy on genus g surface 

The theory on a genus g surface can be analyzed in 
the same way as the theory on torus. The defining 
holonomies of gauge fields are now along 2g elementary 
cycles, i.e. non-contractible loops forming a homology 
basis on the surface. The canonical basis can be made 
by choosing loop pairs {7*,7 2 }, with i-th pair tied to 
i — th "hole" in the surface, i € {1, . . . , <?}, so that only 
intersections of loops are between and 7? (see also 
Ref.Sgi. 

We start with the untwisted sector. As on torus, 



Eq. (C13), the holonomies now lead to labeling of states 
by 



(C26) 



where as in Eq. (C14) the equations of motion lead to 
identifications: 

(C27) 

The I12 invariance to a ■<-> ft allows only the com- 
binations ®i \pi,qi;pi,qi) and ®i{p i rx ,q i a ;p%q 

Pp> QpiPai Qa))- There are (n 2 ) 9 of the former and 

(n 2 ■ n 2 ) 9 — (n 2 ) 9 ]/2 of latter, and since there is only 
one untwisted sector, a total of 



c(") =n 29 '- 

g , untwisted 



,2 S 



1 



(C28) 



Next consider the 2 2s — 1 twisted sectors, labeled by the 
combination of cycles across which there is a twist. As 
discussed in detail in Ref. 22J using Dehn twist transfor- 
mations of the surface, Z2 twists across any combination 
of cycles can be reduced to a twist across a single cycle. 
It is enough therefore to consider a Z2 twist across a sin- 
gle elementary cycle, e.g. the 7* loop. Picking 5-th cycle 
7 1 for the twist means that on an a-flavored copy of the 
surface, E^, the g-th hole is cut open and glued to its 
identical copy E^ along the cut. The resulting surface £ 
has genus 2g — 1. A mirror symmetry of E through the 
g-th hole maps the two flavored copies into each other, 
us : Eg O Eg. We can again introduce a single gauge 
field: 

M^ = iK^ T = { A lf\ (CM) 



(C25) leading to the action 



n 

An 



(bdb + bdb) 



(C30) 
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Since the flat connection B is determined by its 
holonomies, it can be parametrized as 



2g-l 



& = 2^(IH (1) +Xk (2) ) 



2g-l 



(C31) 
(C32) 



using the dual (1-form) basis of the homology basis 



(C33) 



Gauge transformations again lead to X\ , X\ being de- 
fined only up to integers, while the action Eq. ( C30 ) be- 
comes 



25-1 

l = 2™ ]T ( x\x\ 

i=l 



X{x: 



(C34) 



esse ntially 2g— 1 copies of the twisted torus sector theory, 
Eq. (C22). The states are therefore labeled by 



|V> =® l \p\q l ) 1 t = 1...2$-l; 



1 



. . n. 

(C35) 

The action of Z2 here exchanges holes positioned on op- 
posite halves of the surface, i.e. i f-> cr^z, where the 
i = g variables are mapped into themselves. "Diago- 
nal" states have p 1 , q l , i — 1 ... g — 1 quantum numbers 
repeated for i = g + 1 . . . 2g — 1 on the other half of 
the surface, while p 9 ,q 9 are chosen independently; this 
gives (n 2 ) 9 " 1 ■ n 2 — r? 9 invariant states. The "off- 
diagonal" states are symmetrized, and have (p l ,q l ) 7^ 
(p CTE *,f>°' sl ) for i g, while p 9 ,q 9 are independent; this 
gives n 2 ■ (n 2 ) 9 " 1 ^ 2 ) 9 - 1 - l)/2. In total, 



S 



(n) 

g,twisted 



7.2(9" 1) 



1 



(C36) 



Taking into account the number of sectors, Eqs. (C28) 



and (C36 1 give 
n 2g 



c(n) 
^3 



,2» 



± + (2 29 -l)n 29 — ±± (C37) 



2 ' 2 

which for g = 1 reduces to the torus results, 



Eqs. (C15), (C25) 



The final expression for GSD in Eq. (C37) precisely 



matches the one obtained from a simpler physical deriva- 



tion of the melted state, Eq. (26), and it also matches (as 
it must) the GSD result obtained directly from the struc- 
ture of the gauge group G n , Eq. (28). 



Appendix D: GSD of D4 topologically ordered state 
on the three-torus T 3 



To build an exactly solvable model with topological 
order described by discrete (possibly non-Abelian) gauge 



7, Y 



93 



91 

■ ► 




■93 



1 
1 
1 

▲ 



9i 



FIG. 15: A minimal representation of three-torus, having one 
site (dot) of a cubic lattice with periodic boundary conditions 
in all three directions. There are three independent group el- 
ements gi,Q2, <?3 on directed edges (arrows) of this unit-cell, 
the other edges being determined by periodicity. Model with 
topological order with gauge group G in this system contains 
four commuting operators: three force the flux on three cu- 
bic faces to zero, while the fourth acts similarly to a lattice 
gauge transformation on the six directed edges connected to 
the lattice site (dashed lines are periodic images of edges). 



group G in three dimensions with periodic boundary con- 
ditions (i.e. on a three-torus T 3 ), it is enough to consider 
a single unit-cell of a cubic lattice, Fig. [l5j We follow 
the model construction from Refs l49l5T)l There is only 
one lattice site and three directed edges in the unit-cell 
cube, the rest being automatically determined by the pe- 
riodicity. Each edge is assigned a direction (marked by 
arrow in Fig. 15) and a group element, so the degrees of 
freedom are gi G G, i = 1,2,3. The orthonormal basis 
of the Hilbert space is |<?i, 32, <73)- When an edge with 
group element g is considered in direction opposite to its 
arrow, the assigned group element is considered as the 
inverse g~ l . 

The model is defined using operator A (tied to the 
single lattice site) and operators Bij tied to the three 
faces of the unit-cell cube. The B^ operator projects 
the total flux on the ij face of the cube to zero, e.g. 
Bxz enforces the identity ffi • 33 • gf 1 • g^ 1 = 1, where 
direction of edges is taken into account as one goes 
around the XZ square plaquette. (Formally, the zero- 
flux definition leads to the expression Bxz {91,92,93} — 

E 6(hi,gi)6(h2,g 3 )\gi,g2,ga)-) 

h!,h 2 eG 
[hi,h a ]=0 



Next we consider the A g operator, which acts as a 
gauge transformation on the six edges connected to the 
cubic lattice site: Element (/, on edge i is transformed 
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into g ■ gi [p, ■ g^ 1 ] if the edge is directed away [to- 
wards] the site. Note that three of the six edges are peri- 
odic images (dashed lines in Fig. 15). The site operator 



A = l/\G\ J2 g eG therefore acts as 



A\9i,92,9a) = \S919 1 ,99i9 1 ,99z9 l ) ■ (Dl) 



gee 



It is easy to check that A and By operators are all 
projectors and commute with each other. The exactly 
solvable Hamiltonian is therefore given by 



An arbitrary wavefunction 
c(hi,h 2 ,h 3 )\hi,h 2 ,h 3 ), 



is = 
and we consider 

hi,h,2,h 3 £G 

the set of equations for the coefficients c{hi,h 2 ,h 3 ) 
which result from applying Eqs. (D3|. (Note that the 
three By equations just enforce that c(hi,h2,hs) = 
unless hi,h 2l h 3 commute with each other.) 

By solving the equations we find that out of | Z?4 1 3 = 
512 coefficients c(hi,h 2 ,h 3 ) the ground state manifold 
has exactly 92 independent ones, thus the GSD of D4 
topologically ordered state on the three-torus is 92. 



H 



-A-B 



XY 



B 



YZ 



B 



xz ? 



(D2) 



and describes a phase with G topological orderP2HSl xhe 
ground state manifold is given by solving the set of equa- 
tions: 



A\i/>) = \1>) 



(D3) 



Appendix E: Local "surgery" for operators in lattice 
with dislocations 



Here we present the explicit local "surgery" operations 
for redefining local star and plaquette operators, as well 
as other details for the 3d toric code analyzed in the main 



text. All the information is given in Figures 16 to 20 
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FIG. 16: Edge dislocation of an ZX loop, example of edge-1 segment along z. Black dots are removed sites, forming the a/3 
surface. The dislocation line is represented by the outmost array of black sites, and in this edge-1 example does not contain 
star operators. All stabilizers that contain removed spins are repaired by just gluing their neighboring sites together in the 
local neighbor network, and that is how the flavor jump occurs as a string operator which pierces the a/3 plane is constructed 
step-by-step. The exception to the seamless repair are stabilizers that contain spins on the dislocation lines, for which the 
explicit repair "surgery", such that all stabilizers commute, is shown. The two plaquettes (two black triangles) lying in a star 
containing plane orthogonal to the dislocation line are merged into a new plaquette, with a z acting on 9 blue square spins. 
Notice that this stabilizer locally mixes the flavors. The star operator (third black triangle) just loses the removed spin, acting 
on 5 purple diamonds. The surgery is repeated along edge, shifting by a lattice constant (two planes). A new local cubic 
constraint for plaquettes, containing the merged plaquette operator, can be constructed, using er z on the light green spins. Two 
black circles mark the removed local constraints. The corner site of a square-shaped dislocation loop made of edge-1 is on an 
empty site, and therefore needs no special surgery. Locally in total, one plaquette and one constraint are removed, consistent 
with independence of GSD on the length of dislocation lines. 
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FIG. 17: Edge dislocation of an ZX loop, example of edge-2 segment along z. In this case the dislocation line contains star 
operators. See caption of Fig. 16 for details. The two star operators (black triangles) are merged into a new one acting on 13 



purple diamonds. Two plaquettes are replaced by three new ones Bi, B2, -B3, marked by 3 blue, 2 cyan, and 4 purple squares, 
respectively. Three cubic constraints (green circles) are replaced by new ones marked a, b, c, each at a center of a "modified 
cube" formed by neighboring plaquettes. The a, b, c constraints include: two Bi type plaquettes as "cube" top and bottom (in 
total 6 plaquettes); a B3 as a "cube" back side (in total 6); and both B2 and B3 as "cube" sides (in total 7), respectively. 
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FIG. 18: Screw and edge-A dislocation segments of an XY loop (YZ is equivalent). There are no removed sites in this dislocation 
loop, just a shearing of bonds site along y within the loop. All stabilizers are redefined with such a new neighbor network on 
the lattice. 1) The screw dislocation line segment passes through the empty sites of the lattice, while edge-A passes through 
line without star operator sites. An ZX plane plaquette positioned on a screw dislocation line segment (example: rightmost 
black triangle for flavor /?) is repaired by replacement by a 9 site a z operator B s \ (cyan squares). The cubic constraints on 
this segment (rightmost black circle for flavor a) now use the repaired plaquette operators on front and back cube sides. 2) For 
each spin site R on the edge-A segment, two plaquettes (black triangles) are merged into one new B„2 (blue squares, mixing 
flavors), while the star positioned below R is edited by removing its action on spin at R. Two local cubic constraints on either 
side of edge-A sites are merged into a single constraint containing two B S 2 plaquettes. Green circles mark plaquettes which are 
involved, together with B a i and B S 2, in the cubic constraint (black circle) at the loop corner. The corner should be included 
in the screw segment. In total locally, the number of stabilizers and constraints is unchanged, consistent with independence of 
GSD on the size of dislocation loops. 
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FIG. 19: Edge-B dislocation segment of an XY loop (YZ is equivalent), lying along x (inequivalent to edge-A due to b — y). 
See caption of Fig. |18| The edge-B does not pass through star operator sites. The plaquettes labeled by spins on the dislocation 
line sites (example: black triangle for flavor /3) have to lose one spin operator (blue squares). The plaquettes on empty sites of 
edge-B (examples are black circles) are replaced by two (purple and cyan squares). Pairs of star operators below edge-B (black 
triangles) are merged (green diamonds). The corner (marked with a) should be included into the edge-B segment. In total 
locally, the number of stabilizers (stars and plaquettes) and local cubic constraints is unchanged, consistent with independence 
of GSD on the size of dislocation loops. The cubic constraint belonging to corner (a) includes plaquettes labeled by green circles. 
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FIG. 20: Dislocation loop linking, example of ZX and XY loops. After the screw loop is created by editing the local neighbor 
network, the sites for the ZX loop should be removed, and network edited again per rules of Fig. [16] A single plaquette in 
the screw plane (bla ck triangle in circle) must be completely removed. The cubic constraint at that same site now contains B\ 
plaquettes (Fig. 16 1 and the ones marked by green circles. 
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